Draft version May 2 1 , 2008 

Preprint typeset using I4TgX style emulateapj v. 08/13/06 



A MULTIPLICITY CENSUS OF YOUNG STARS IN CHAMAELEON I 

David Lafreniere 1 , Ray Jayawardhana 1 , Alexis Brandeker 2 , Mirza Ahmic 1 , Marten H. van Kerkwijk 1 

Draft version May 21, 2008 



OO 
O 
O 

(N 



Si 

6 

(N 
> 

in 
p 

cn 
o 

OO 

o 



X 



ABSTRACT 

We present the results of a multiplicity survey of 126 stars spanning ~0.1— 3 M Q in the ^2-Myr-old 
Chamaeleon I star-forming region, based on adaptive optics imaging with the ESO Very Large Telescope. 
Our observations have revealed 30 binaries and 6 triples, of which 19 and 4, respectively, are new discoveries. 
The overall multiplicity fraction we find for Cha I (~30%) is similar to those reported for other dispersed young 
associations, but significantly higher than seen in denser clusters and the field, for comparable samples. Both 
the frequency and the maximum separation of Cha I binaries decline with decreasing mass, while the mass 
ratios approach unity; conversely, tighter pairs are more likely to be equal mass. We confirm that brown dwarf 
companions to stars are rare, even at young ages at wide separations. Based on follow-up spectroscopy of two 
low-mass substellar companion candidates, we conclude that both are likely background stars. The overall 
multiplicity fraction in Cha I is in rough agreement with numerical simulations of cloud collapse and fragmen- 
tation, but its observed mass dependence is less steep than predicted. The paucity of higher-order multiples, in 
particular, provides a stringent constraint on the simulations, and seems to indicate a low level of turbulence in 
the prestellar cores in Cha I. 

Subject headings: Stars: formation — binaries: general — stars: low-mass, brown dwarfs 



1. INTRODUCTION 

The ubiquity of multiple stellar systems, both in the field 
(e.g. Duquennoy & Mayor 1991) and in star-forming regions 
(e.g. Ghez et al. 1997), clearly indicates that multiplicity is a 
principal channel of star formation. Accordingly, the multi- 
plicity properties of a stellar population - such as the fraction 
of binary or higher-order multiple systems, the orbital separa- 
tion and mass ratio distributions of companions, and the de- 
pendence of these on stellar mass - should be good tracers of 
the specific mechanism by which the stars formed, and possi- 
bly also of the characteristics of their birth environment (e.g. 
Goodwin et al. 2004a; Bate et al. 2003). Beyond this primor- 
dial imprint, many of these properties can be further modified 
by dynamical interactions or decay of the newborn systems 
over a few million years (e.g. Sterzik & Durisen 1998, 2003; 
Kroupa 1995, 1998), a process that may also be environment 
dependent and that could introduce spatial variations of mul- 
tiplicity properties within a given cluster. Studies of stellar 
multiplicity can thus provide valuable insight into the global 
process of star formation and early evolution, especially if 
multiplicity properties can be compared across a wide range 
of stellar mass, for different ages, and for different environ- 
ments. 

The multiplicity properties of stellar systems in the Solar 
neighborhood have been studied extensively (e.g. Duquen- 
noy & Mayor 1991; Fischer & Marcy 1992; Delfosse et al. 
2004), at least for ^0.1-2 M stars. The picture is not yet as 
complete for substellar objects, but progress is being made 
rapidly (see Burgasser et al. 2007, and references therein). 
For pre-main-sequence stars, since they usually lie more than 
100 pc away from Earth, observational advances have been 
closely tied to the development of large aperture telescopes 
and high angular resolution imaging techniques. The first 
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speckle imaging surveys of a few star-forming regions, carried 
out in the early 1990s, revealed binary frequencies that ap- 
peared to be significantly higher than in the field (Ghez et al. 
1993; Leinert et al. 1993). This early realization prompted 
observations of more young clusters, and the ensuing results 
suggested that multiplicity may vary with the environment 
(Ghez et al. 1997; Padgett et al. 1997; Petr et al. 1998; Kohler 
& Leinert 1998; Simon et al. 1999; Duchene et al. 1999; 
Kohler et al. 2000). While dispersed Taurus-like associations 
show an overabundance of binaries by a factor ^2 in compar- 
ison to the field, the binary fraction in dense Trapezium-like 
clusters appears to be in good agreement with the field value 
(Petr et al. 1998; Simon et al. 1999). Whether this dichotomy 
is a direct outcome of the star formation process or the result 
of stronger dynamical evolution in denser environments re- 
mains unclear (Petr et al. 1998; Duchene et al. 1999; Mathieu 
et al. 2000; Brandeker et al. 2006). 

Typically, owing to the brightness requirements inherent 
to speckle imaging, many of the early surveys neglected the 
fainter, very low mass (VLM) members of young associations 
and clusters. However, a number of recent multiplicity stud- 
ies have specifically targeted brown dwarfs and VLM stars in 
star-forming regions (e.g. Rraus et al. 2005, 2006; Bouy et al. 
2006; Ahmic et al. 2007; Konopacky et al. 2007). The young 
VLM binaries found in these surveys share several character- 
istics with their counterparts in the field: lower binary frac- 
tion than for higher mass stars, lower average orbital separa- 
tion, and a preference for mass ratios near unity (Close et al. 
2003; Burgasser et al. 2003; Delfosse et al. 2004; Burgasser 
et al. 2007). But, the current observations also hint at a higher 
incidence of wide (>30 AU) VLM binaries in star-forming 
regions in comparison with the field (e.g. Konopacky et al. 
2007; Jayawardhana & Ivanov 2006). 

The compendium of high angular resolution imaging data 
of pre-main-sequence stars now covers a wide range of en- 
vironments, ages and masses. However, the various surveys 
often have different resolution and contrast limits, and their 
samples are drawn from different mass regimes. As a result, 
determining the mass dependence of multiplicity and compar- 
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ing the properties of different star-forming environments are 
no easy tasks. This is rather unfortunate, especially given that 
recent theoretical efforts, which have focused on modeling the 
collapse and fragmentation of prestellar cores and the early 
dynamical interactions of small-N stellar clusters (e.g. Sterzik 
& Durisen 2003; Delgado-Donate et al. 2004; Goodwin et al. 
2004a), have produced results that could be compared directly 
with observable quantities. Thus, we have embarked upon 
an ambitious program to obtain a comprehensive and detailed 
picture of the multiplicity properties of pre-mains-sequence 
stars across a wide mass range, based on a dataset of homoge- 
neous quality and sensitivity for several star-forming regions. 

Here we present the results of an adaptive optics (AO) 
multiplicity survey of 126 stars spanning ~0. 1—3 M Q in the 
nearby (160 pc; Whittet et al. 1997), young (2 Myr; Luhman 
2004 and references therein) Chamaeleon I (hereafter Cha 
I) cloud. Ahmic et al. (2007) presented an analysis of the 
subsample of VLM targets from this survey while Damjanov 
et al. (2007) investigated circumstellar disks of Cha I mem- 
bers, including the effect of companions on disks. The present 
paper reports the results for the entire sample observed, with 
a particular emphasis on the multiplicity properties as a func- 
tion of stellar mass. 

The structure of the paper is as follows. The target sample 
and the observations are described in §2, and the data analy- 
sis performed is detailed in §3. The multiple systems identi- 
fied and their measured properties are presented in §4. The 
multiplicity characteristics of our sample, including the mul- 
tiplicity fraction and its dependence on mass, the mass ratio 
of companions as a function of primary mass and orbital sep- 
aration, the binary separations as a function of mass, and the 
stability of the triple systems, are discussed in §5 in relation to 
other stellar populations and to numerical simulations of star 
formation. Finally, concluding remarks follow in §6. 

2. OBSERVATIONS AND DATA REDUCTION 

The target sample was built upon the recent compilation of 
Cha I members made by Luhman (2004). The membership of 
these stars was established based on space velocity, H-R di- 
agram position, infrared excess emission, lithium abundance, 
and stellar activity indicators. The only selection criterion ap- 
plied to the list of Luhman (2004) was a magnitude cut at 
K=13, the faint magnitude limit of the IR wave sensor of the 
AO system used for the observations; this removed only 8 out 
of the 158 known Cha I members and does not introduce a 
significant bias. We have observed 126 out of the remaining 
150 Cha I members, our sample is thus a good representation 
of the entire population of Cha I. The basic properties of the 
stars observed are indicated in Table l. 3 The TT-band magni- 
tudes, spectral types and effective temperatures of the targets 
were taken from Luhman (2004). 

The observations were obtained with the Nasmyth Adap- 
tive Optics System (NAOS) and the CONICA camera at the 
8.2m telescope Yepun, part of the European Southern Obser- 
vatory's Very Large Telescope (program 076.C-0579). The 
data were acquired in visitor mode over three consecutive 
nights starting on March 24, 2006, with the exception of the 



source Cha Ha 2, which was observed in service mode on 
2005 March 25 as part of a preliminary run (program 075. C- 
0042). The high-resolution lens of CONICA was used for 
all observations, yielding a pixel scale of 13.27 mas and a 
field of view of 13.6" x 13.6". For the majority of targets, 
the N90C10 dichroic was used to direct 90% of the light to 
the AO wave front sensor and 10% to the science camera, 
but whenever possible the VIS dichroic was used to take ad- 
vantage of its higher throughput toward CONICA at near-IR 
wavelengths. All targets have been observed in the K band us- 
ing either the K s or the narrow-band 2.17 fim filter, and most 
targets have been further observed in the H band using either 
the H or the narrow-band 1.64 /im filter. A typical imaging 
sequence consisted of six exposures divided over two dither 
positions separated by ~5.5". Observing parameters for each 
target are indicated in Table 1. The data were reduced in 
a standard manner using the NAOS-CONICA pipeline. As 
an indication of image quality, the TT-band point-spread func- 
tion (PSF) full-width-at-half-maximum (FWHM) is indicated 
in Table 1 . The image quality was highly variable over the 
three nights of observations, and the median TT-band FWHM 
observed is 0.086". For reference, the /f-band diffraction- 
limited FWHM is 0.056". 

The images of the stars CHXR 18N and T54 showed 
faint candidate companions that, if truly bound, would have 
masses well below the substellar limit. However, being faint 
and located at relatively large separations from their primary 
(2.7"and 2.4", resp.), the probabilities that these sources are 
actually unrelated background stars are significant. Based on 
the analysis presented in §4.2, these probabilities are 0.33 
and 0.40, respectively. To assess whether or not these can- 
didates are truly bound companions and to clarify their na- 
ture, we have obtained follow-up imaging and spectroscopy 
at the VLT on 2007 May 13 for T54 [278.C-5070(A)] and on 
2007 June 18 for CHXR 18N [279.C-5017(A)]. The imaging 
follow-up was obtained in 7, H, and K s using the same instru- 
mental configuration as our main survey and data reduction 
was carried out as described above. For T54, the total expo- 
sure times were 7.5 min, 9 min, and 4.8 min in J, H, and K s , 
respectively. For CHXR 18N, they were 9 min, 2.8 min, and 
2 min, respectively. The spectroscopic observations were ac- 
quired using the VIS wavefront sensor, a 40" x 0.1 72" slit, 
and the S27-4-SH spectroscopic mode of NACO, covering 
the wavelength range 1.37-1.84 /im. Exposures at two dif- 
ferent positions along the slit were obtained two times in an 
ABBA pattern, with each position observed 63 x 8 s for T54 
and 9 x 1 8 s for CHXR 1 8N, totaling 4 x 63 x 8 s = 34 min in- 
tegration on T54 and 4x9x18 s = 11 min on CHXR 18N. 
For the binary T54 and its candidate companion as well as 
for CHXR 18N and its candidate companion, the slit was 
aligned with the primary and candidate companion; for T54, 
this means that the tight secondary component was out of the 
slit. Immediately after each observation, a calibration star was 
observed with the same settings and at a similar airmass to 
correct for atmospheric and instrumental transmission. The 
B4 star HIP 53024 was used for T54 and the solar-like G3 
star HIP 88299 was used for CHXR 18N. 



3 Six other stars were observed during the observing run but are not con- thesg stars ^ omitted from TaWe , 

sidered in this paper as the AO correction achieved was extremely poor 
(Ced 110 IRS 6, Cha Ha 10, Cha Ha 11, IRN, ISO 217, and ISO 220); 
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TABLE 1 

Target properties and observation log 





RA 


DEC 


Ki 


Spectral 


r eff a 


Mass^ 1 




(s) Dich. c 


FWHM d 


Detection limits (K s ) 


Name 


(J2000.0) 


(J2000.0) 


\ LLLa 5/ 


type 


(K) 




K s 


H 


("1 


0.1" 


0.25" 


0.5" 


1.0" 


CHSM 17173 


11 10 22.27 


-76 25 13.8 


12.45 


M8 


2710 


0.03 


240 


1 


0.131 


13.7 


15.2 


15.8 


15.8 


Cha Ha 1 


11 07 16.69 


-77 35 53.3 


12.17 


M7.75 


2752 


0.06 


120 


120 1 


0.144 


13.2 


14.2 


14.5 


14.5 


Cha Ha 7 


11 07 37.76 


-77 35 30.8 


12.42 


M7.75 


2752 


0.06 


240 


1 


0.133 


13.5 


14.9 


15.5 


15.6 


Cha Ha 12 


1 1 06 38.00 


-77 43 09.1 


11.81 


M6.5 


2935 


0.12 


120 


120 1 


0.119 


13.1 


14.5 


15.0 


15.0 


ISO 138 


11 08 18.50 


-77 30 40.8 


13.04 


M6.5 


2935 


0.12 


240 


1 


0.121 


14.2 


15.0 


15.2 


15.1 


CHSM 1982 


11 04 10.60 


-76 12 49.0 


12.12 


M6 


2990 


0.14 


120 


120 1 


0.092 


13.9 


15.0 


15.2 


15.2 


ISO 252 


11 10 41.41 


-77 20 48.0 


12.27 


M6 


2990 


0.14 


120 


1 


0.135 


13.4 


14.7 


15.1 


15.1 


CHSM 10862 


11 07 46.56 


-76 15 17.5 


12.33 


M5.75 


3024 


0.14 


120 


1 


0.123 


13.6 


14.6 


15.0 


15.0 


Cha Ha 6 


11 08 39.52 


-77 34 16.7 


11.04 


M5.75 


3024 


0.14 


120 


120 1 


0.105 


12.9 


14.9 


15.5 


15.5 


Hn 13 


11 10 55.97 


-76 45 32.6 


9.91 


M5.75 


3024 


0.14 


120 


108 1 


0.072 


12.3 


14.3 


15.3 


15.9 


ISO 147 


11 08 26.51 


-77 15 55.1 


12.35 


M5.75 


3024 


0.14 


240 


1 


0.144 


13.4 


14.6 


15.2 


15.3 


CHXR 84 


11 12 03.27 


-76 37 03.4 


10.78 


M5.5 


3058 


0.15 


120 


120 1 


0.091 


13.0 


15.2 


16.0 


16.1 


Cha Ha 13 


11 08 17.03 


-77 44 11.8 


10.67 


M5.5 


3058 


0.15 


120 


120 1 


0.096 


12.9 


15.3 


16.2 


16.3 


Cha Ha 3 


11 07 52.26 


-77 36 57.0 


11.10 


M5.5 


3058 


0.15 


120 


120 1 


0.096 


13.3 


15.7 


16.4 


16.4 


Cha Ha 4 


11 08 18.96 


-77 39 17.0 


11.02 


M5.5 


3058 


0.15 


120 


120 1 


0.086 


13.3 


15.5 


16.1 


16.1 


Cha Ha 5 


11 08 24.11 


-77 41 47.4 


10.71 


M5.5 


3058 


0.15 


120 


120 1 


0.113 


12.5 


14.7 


15.6 


15.5 


Cha Ha 9 


11 07 18.61 


-77 32 51.7 


11.80 


M5.5 


3058 


0.15 


120 


120 1 


0.099 


13.4 


14.9 


15.4 


15.4 


Hn 12W 


11 10 28.52 


-77 16 59.6 


10.78 


M5.5 


3058 


0.15 


120 


120 1 


0.095 


12.7 


14.7 


15.4 


15.5 


ISO 165 


1 1 08 54.97 


-76 32 41.1 


11.44 


M5.5 


3058 


0.15 


120 


120 1 


0.132 


12.5 


13.7 


14.4 


14.5 


ISO 235 


11 10 07.85 


-77 27 48.1 


11.34 


M5.5 


3058 


0.15 


120 


120 1 


0.103 


13.0 


14.8 


15.7 


15.8 


ISO 28 


11 03 41.87 


-77 26 52.0 


11.69 


M5.5 


3058 


0.15 


120 


120 1 


0.098 


13.3 


14.8 


15.3 


15.3 


CHXR 15 


11 05 43.00 


-77 26 51.8 


10.23 


M5.25 


3091 


0.16 


30 


30 1 


0.086 


12.6 


14.2 


14.5 


14.5 


CHXR 78C 


1 1 08 54.22 


-77 32 11.6 


11.22 


M5.25 


3091 


0.16 


30 


30 1 


0.112 


12.9 


13.6 


13.5 


13.5 


Cha Ha 2<= 


1 1 07 42.45 


-77 33 59.4 


10.68 


M5.25 


3091 


0.16 


1080 


1 


0.087 


12.4 


13.9 


14.5 


14.9 


C7-1 


1 1 09 42.60 


-77 25 57.9 


10.55 


M5 


3125 


0.17 


30 


30 1 


0.126 


12.0 


13.5 


13.7 


13.6 


Hn2 


11 03 47.64 


-77 19 56.3 


9.99 


M5 


3125 


0.17 


30 


30 1 


0.072 


12.7 


14.9 


15.5 


15.5 


ISO 143 


11 08 22.38 


-77 30 27.7 


11.10 


M5 


3125 


0.17 


30 


30 1 


0.098 


12.9 


14.1 


14.3 


14.3 


T50 


11 12 09.85 


-76 34 36.6 


9.84 


M5 


3125 


0.17 


30 


30 1 


0.069 


ion 
12. V 


1 j. j 


1 J. 9 


1 J. 9 


Hn7 


11 09 05.13 


-77 09 58.1 


10.96 


M4.75 


3161 


0.17 


30 


30 1 


0.090 


13. U 


1 A A 

14.4 


1 A & 

14. 


1 A Z 

14. j 


ISO 250 


1 1 10 36.45 


-77 22 13.2 


10.67 


M4.75 


3161 


0.17 


30 


30 1 


0.100 


IZ.u 


1 j.y 


1 A 1 

14.1 


1 A 1 

14.1 


B35 


11 07 21.43 


-77 22 11.8 


10.93 


- 


- 




30 


1 


0.113 


12.5 


13.6 


13.8 


13.8 


Hn5 


11 06 41.81 


-76 35 49.0 


10.12 


M4.5 


3198 


0.18 


30 


30 1 


0.081 


12.4 


14.3 


14.9 


15.0 


ISO 274 


1111 22.61 


-77 05 53.9 


10.69 


M4.5 


3198 


0.18 


30 


30 1 


0.092 


12.7 


14.4 


14.8 


14.8 


T12 


11 02 55.05 


-77 21 50.8 


10.45 


M4.5 


3198 


0.18 


30 


30 1 


0.089 


12.8 


14.8 


15.2 


15.1 


T55 


11 13 33.57 


-76 35 37.4 


10.73 


M4.5 


3198 


0.18 


30 


30 1 


0.098 


13.0 


14.7 


14.8 


14.7 


CHXR 60 


11 13 29.71 


-76 29 01.2 


10.58 


M4.25 


3234 


0.19 


30 


30 1 


0.091 


12.4 


13.9 


14.2 


14.1 


CHXR 74 


11 06 57.33 


-77 42 10.7 


10.21 


M4.25 


3234 


0.19 


30 


30 1 


0.095 


11.7 


13.0 


13.6 


13.7 


CHXR 76 


11 07 35.19 


-77 34 49.3 


10.95 


M4.25 


3234 


0.19 


30 


30 1 


0.094 


13.0 


14.5 


14.8 


14.7 


2M1110-7722 h 


11 10 34.81 


-77 22 05.3 


10.03 


M4 


3270 


0.20 


30 


30 1 


0.083 


12.2 


14.3 


15.1 


15.1 


Hn21W 


11 14 24.54 


-77 33 06.2 


10.65 


M4 


3270 


0.20 


120 


120 1 


0.091 


12.5 


14.7 


15.8 


16.0 


ISO 52 


11 04 42.58 


-77 41 57.1 


10.64 


M4 


3270 


0.20 


30 


30 1 


0.097 


12.6 


14.1 


14.5 


14.4 


CHXR 62 


11 14 15.65 


-76 27 36.4 


10.12 


M3.75 


3306 


0.22 


30 


30 1 


0.087 


11.9 


13.2 


13.6 


13.7 


T10 


1 1 00 40.22 


-76 19 28.1 


10.87 


M3.75 


3306 


0.22 


30 


30 1 


0.094 


13.1 


14.3 


14.4 


14.5 


T34 


11 08 16.49 


-77 44 37.2 


10.02 


M3.75 


3306 


0.22 


30 


30 1 


0.072 


12.7 


14.6 


15.0 


15.1 


CHXR 12 


11 03 56.83 


-77 21 33.0 


9.71 


M3.5 


3342 


0.23 


30 


30 1 


0.082 


12.1 


14.1 


14.8 


14.9 


CHXR 22E 


11 07 13.30 


-77 43 49.9 


9.99 


M3.5 


3342 


0.23 


30 


30 1 


0.098 


11.7 


13.9 


14.7 


14.6 


CHXR 26 e 


11 07 36.87 


-77 33 33.5 


9.35 


M3.5 


3342 


0.23 


30 


30 1 


0.085 


11.5 


14.0 


14.8 


14.7 


Hn 18 


11 13 24.46 


-76 29 22.7 


10.80 


M3.5 


3342 


0.23 


30 


30 1 


0.081 


13.4 


15.2 


15.5 


15.4 


B43 


11 09 47.42 


-77 26 29.1 


10.24 


M3.25 


3379 


0.24 


30 


30 1 


0.090 


12.2 


14.3 


15.1 


15.1 


CHXR 73 k 


11 06 28.77 


-77 37 33.2 


10.70 


M3.25 


3379 


0.24 


30 


30 1 


0.099 


12.6 


14.1 


14.5 


14.5 


Hn 10E 


11 09 46.21 


-76 34 46.4 


10.05 


M3.25 


3379 


0.24 


30 


30 1 


0.076 


12.6 


14.5 


15.0 


15.1 


Hn4 


11 05 14.67 


-77 11 29.1 


9.61 


M3.25 


3379 


0.24 


30 


30 1 


0.081 


12.1 


13.9 


14.3 


14.5 


T5 


10 57 42.20 


-76 59 35.7 


9.25 


M3.25 


3379 


0.24 


30 


30 1 


0.072 


12.2 


13.7 


14.3 


14.9 


CHXR 21 


11 07 11.49 


-77 46 39.4 


9.66 


M3 


3415 


0.26 


30 


30 1 


0.098 


11.6 


13.8 


14.6 


14.7 


CHXR 71 


11 02 32.65 


-77 29 13.0 


10.13 


M3 


3415 


0.26 


30 


30 1 


0.070 


12.9 


14.8 


15.1 


15.1 


T16 


11 04 57.01 


-77 15 56.9 


10.41 


M3 


3415 


0.26 


30 


30 1 


0.086 


12.7 


14.6 


15.2 


15.1 


T22 


11 06 43.47 


-77 26 34.4 


9.39 


M3 


3415 


0.26 


30 


30 1 


0.076 


12.0 


14.5 


15.4 


15.4 


B53 


11 14 50.32 


-77 33 39.0 


9.55 


M2.75 


3451 


0.27 


30 


30 1 


0.079 


11.9 


14.0 


14.7 


14.8 


CHXR 57 


11 13 20.13 


-77 01 04.5 


10.01 


M2.75 


3451 


0.27 


30 


30 1 


0.069 


13.2 


15.5 


15.9 


15.9 


CHXR 59 


11 13 27.37 


-76 34 16.6 


9.63 


M2.75 


3451 


0.27 


30 


30 1 


0.078 


12.1 


13.7 


14.3 


14.5 


Cam2-19 


11 06 15.45 


-77 37 50.1 


10.25 


M2.75 


3451 


0.27 


30 


30 1 


0.094 


12.3 


14.2 


14.8 


14.6 


CHXR 48 


11 11 34.75 


-76 36 21.1 


9.80 


M2.5 


3488 


0.29 


30 


30 1 


0.078 


12.4 


14.6 


15.2 


15.2 


T25 


11 07 19.15 


-76 03 04.8 


9.77 


M2.5 


3488 


0.29 


30 


30 1 


0.075 


12.5 


14.9 


15.5 


15.5 


T30 


11 07 58.09 


-77 42 41.3 


9.89 


M2.5 


3488 


0.29 


30 


30 1 


0.079 


12.4 


14.6 


15.2 


15.2 


CHXR 9C 


11 01 18.75 


-76 27 02.5 


8.99 


M2.25 


3524 


0.32 


30 


30 1 


0.082 


11.4 


13.3 


13.5 


13.5 


T39S 


11 09 11.72 


-77 29 12.5 


8.96 


M2 


3560 


0.35 


48 


30 1 


0.079 


11.4 


13.9 


14.6 


14.5 


T43 


11 09 54.08 


-76 29 25.3 


9.25 


M2 


3560 


0.35 


30 


30 1 


0.072 


12.2 


14.7 


15.4 


15.4 


T47 


11 10 49.60 


-77 17 51.8 


9.18 


M2 


3560 


0.35 


30 


30 1 


0.072 


11.6 


13.9 


14.9 


14.9 


T49 


1111 39.66 


-76 20 15.2 


8.87 


M2 


3560 


0.35 


48 


30 1 


0.090 


10.9 


13.1 


13.9 


14.0 


CHXR 14S 


11 04 52.85 


-76 25 51.5 


9.75 


M1.75 


3596 


0.37 


30 


30 1 


0.078 


12.1 


14.4 


15.1 


15.0 


T20 


11 05 52.61 


-76 18 25.6 


9.34 


M1.5 


3632 


0.40 


30 


30 1 


0.069 


12.3 


14.8 


15.7 


15.7 


T23 


11 06 59.07 


-77 18 53.6 


10.00 


M1.5 


3632 


0.40 


30 


30 1 


0.070 


12.9 


15.1 


15.6 


15.5 


Cl-6 


1 1 09 22.67 


-76 34 32.0 


8.67 


M1.25 


3669 


0.42 


30 


30 1 


0.083 


10.7 


13.3 


14.5 


14.6 


CHXR 30B 


11 07 57.31 


-77 17 26.2 


9.95 


M1.25 


3669 


0.42 


30 


30 1 


0.076 


12.3 


14.6 


15.5 


15.5 
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TABLE 1 — Continued 





RA 


DEC 




Spectral 


7* e ff a 


Mass 11 


t 

K s 


(s) 


Dich. c 


FWHM d 


Detection limits (K s ) 


Name 


(J2000.0) 


(J2000.0) 


(mag) 




(K) 




M 




(") 


0.1" 


0.25" 


0.5" 


1.0" 


CHXR 40 


1 1 09 40.07 


-76 28 39.2 


8.96 


M1.25 


3669 


0.42 


48 


30 


1 


0.090 


10.9 


12.7 


13.4 


13.6 


CHXR 79 


11 09 18.13 


-76 30 29.2 


9.06 


M1.25 


3669 


0.42 


30 


30 


1 


0.086 


11.3 


13.6 


14.6 


14.7 


ISO 126 


11 08 02.98 


-77 38 42.6 


8.30 


M1.25 


3669 


0.42 


48 


30 


1 


0.079 


10.6 


12.8 


13.9 


14.4 


CHXR 54 


11 12 42.10 


-76 58 40.0 


9.50 


Ml 


3705 


0.44 


30 


30 


1 


0.075 


12.2 


14.7 


1 £ A 
19.4 


1 £ A 

15.4 


T27 e 


11 07 28.26 


-76 52 11.9 


9.52 


Ml 


3705 


0.44 


30 


30 


1 


0.079 


11.9 


13.8 


14.4 


14.6 


T28 


11 07 43.66 


-77 39 41.1 


8.26 


Ml 


3705 


0.44 


48 


30 


1 


0.072 


10.7 


13.1 


1 A 1 

14.3 


1/10 

14.0 


T48 


11 10 53.34 


-76 34 32.0 


10.04 


Ml 


3705 


0.44 


30 


30 


1 


0.078 


1Z.0 


1 A ^ 

14. J 


1C I 

19.1 


1 ^ 1 
19.1 


T53 


11 12 30.93 


-76 44 24.1 


9.12 


Ml 


3705 


0.44 


30 


30 


1 


0.069 


1 9 ft 
1Z.U 


1 A 

14.0 


19. / 


1 ^ c 
19.5 


T24 


11 07 12.07 


-76 32 23.2 


9.38 


M0.5 


3778 


0.48 


30 


30 


1 


0.076 


1 9 ft 
iz.u 


1 A 

14.0 


1 ^ A 
19.4 


1 ^ ^ 

19.9 


T3 C 


10 55 59.73 


-77 24 39.9 


8.69 


M0.5 


3778 


0.48 


30 


30 


1 


0.106 


1 n 9 
1U.Z 


12. 5 


1 A 1 

14.1 


1 A 6. 

14.0 


T38 


11 08 54.64 


-77 02 13.0 


9.46 


M0.5 


3778 


0.48 


30 


- 


1 


0.171 


i ft 9. 
1U. J 


1 1 A 

1 1.4 


1 9 ft 
IZ.U 


19 1 

lz.1 


T4 


10 56 30.45 


-77 11 39.3 


8.63 


M0.5 


3778 


0.48 


30 


- 


1 


0.075 


1 i n 
1 l.U 


1 9 n 
1 J. / 


Ij.U 


1 <C 9 
Ij.Z 


T56 


11 17 37.01 


-77 04 38.1 


9.23 


M0.5 


3778 


0.48 


30 


30 


1 


0.069 


19 9 
IZ.Z 


1 A 

14.0 


1 ^ A 
lj.4 


1 ^ A 
lj.4 


CI -25 


11 09 41.93 


-76 34 58.4 


10.00 


- 


- 


0.49 f 


30 


30 


1 


0.136 


11.2 


12.6 


13.1 


13.1 


CHXR 33 


11 08 40.69 


-76 36 07.8 


9.28 


MO 


3850 


0.51 


30 


30 


1 


0.069 


12.1 


14.6 


15.7 


15.7 


CHXR 49NE 


1111 54.00 


-76 19 31.1 


9.23 


MO 


3850 


0.51 


30 


30 


1 


0.069 


12.3 


14.2 


14.7 


14.9 


T45A 


11 10 04.69 


-76 35 45.3 


9.24 


MO 


3850 


0.51 


30 


30 


1 


0.069 


11.9 


14.5 


15.7 


15.9 


CHXR 14N 


1104 51.00 


-76 25 24.1 


9.60 


K8 


3955 


0.55 


30 


30 


1 


0.079 


12.0 


14.4 


15.2 


15.2 


CHXR 30A 


1108 00.03 


-77 17 30.5 


9.09 


K8 


3955 


0.55 


30 


30 


1 


0.085 


1 1 Q 

1 1.3 


13.5 


14.3 


14.4 


CHXR 68A 


11 18 20.24 


-76 21 57.6 


8.87 


K8 


3955 


0.55 


48 


30 


1 


0.076 


11.4 


13.4 


14.1 


14.3 


Hn 11 


11 10 03.69 


-76 33 29.2 


9.44 


K8 


3955 


0.55 


30 


30 


1 


0.085 


1 1.7 


14.1 


15.1 


15.1 


T318 


11 08 01.49 


-77 42 28.9 


6.96 


K8 


3955 


0.55 


30 


30 


1 


0.075 


9.8 


12.0 


13.4 


14.6 


T35 


11 08 39.05 


-77 16 04.2 


9.11 


K8 


3955 


0.55 


30 


30 


1 


0.069 


11.8 


14.1 


15.1 


15.3 


T45 e 


11 09 58.74 


-77 37 08.9 


7.97 


K8 


3955 


0.55 


30 


30 


1 


0.089 


9.9 


12.0 


13.3 


14.0 


T46 


1 1 10 07.04 


-76 29 37.7 


8.45 


K8 


3955 


0.55 


48 


30 


1 


0.069 


11. z 


IT 11 

19.0 


1 ^ ft 


19. / 


T7 


10 59 01.09 


-77 22 40.7 


8.62 


K8 


3955 


0.55 


30 


30 


1 


0.072 


11.2 


13.9 


15.1 


15.3 


CHXR 37 


11 09 17.70 


-76 27 57.8 


8.70 


K7 


4060 


0.59 


30 


30 


1 


0.070 


10.8 


13.9 


15.3 


19.0 


Cam2-42 


11 09 37.78 


-77 10 41.1 


9.16 


K7 


4060 


0.59 


30 


30 


1 


0.090 


11. j 


1 9. 1 


14. o 


11 o 
14. y 


Cl-2 


11 09 55.06 


-76 32 41.0 


9.67 


- 


- 


0.64 f 


30 


30 


1 


0.108 


11.1 


12.9 


13.8 


13.8 


CHXR 18N 


11 11 46.32 


-76 20 09.2 


7.77 


K6 


4205 


0.65 


30 


- 


1 


0.072 


10.3 


12.7 


14.3 


15.3 


CHXR 20 


1106 45.10 


-77 27 02.3 


8.88 


K6 


4205 


0.65 


30 


30 


1 


0.095 


10.7 


12.7 


13.6 


13.7 


CHXR 28<= 


11 07 55.89 


-77 27 25.8 


7.69 


K6 


4205 


0.65 


30 


30 


1 


0.092 


9.2 


11.0 


12.5 


13.6 


Til 


11 02 24.91 


-77 33 35.7 


8.20 


K6 


4205 


0.65 


48 


30 


1 


0.082 


10.7 


13.6 


15.2 


15.7 


T40 


1109 23.79 


-76 23 20.8 


8.24 


K6 


4205 


0.65 


48 


30 


1 


0.069 


11.0 


13.7 


15.3 


16.0 


ISO 237 


11 10 11.42 


-76 35 29.3 


8.62 


K5.5 


4278 


0.68 


30 


30 


1 


0.079 


1 1 A 
1 l.U 


ii t 
13.6 


1 a n 

14. y 


1 C A 

15.4 


T14 


1104 09.09 


-76 27 19.4 


8.66 


K5 


4350 


0.71 


30 


30 


1 


0.069 


1 1.6 


14.3 


15.5 


15.8 


T42 


11 09 53.41 


-76 34 25.5 


6.46 


K5 


4350 


0.71 


30 


30 


1 


0.107 


1 Q 


Q S 
y.o 


I 1 9 

I I .Z 


1 1 Q 

ii.y 


T44 


11 10 00.11 


-76 34 57.9 


6.08 


K5 


4350 


0.71 


48 


- 


1 


0.121 


7 2 


Q 9 
y.z 


10 7 


115 


CHXR 55 


11 12 43.00 


-76 37 04.9 


9.29 


K4.5 


4470 


0.77 


30 


30 


1 


0.069 


1 9 ^ 
1Z. j 


1 S 1 
19.1 


i f\ ft 


i f\ ft 

1U.U 


T51 e 


11 12 24.41 


-76 37 06.4 


8.00 


K3.5 


4660 


0.88 


48 


30 


1 


0.092 


Q 7 
y. 1 


116 


13 3 


14 5 


CHXR 47<= 


11 10 38.02 


-77 32 39.9 


8.28 


K3 


4730 


0.96 


48 


30 


1 


0.091 


10.5 


12 4 


13 5 


14 3 


T8 


10 59 06.99 


-77 01 40.4 


7.31 


K2 


4900 


1.30 


31 


- 




0.078 


9.8 


12.0 


13.4 


14.9 


T6 e 


10 58 16.77 


-77 17 17.1 


7.76 


K0 


5250 


1.96 


30 


30 


1 


0.072 


10.0 


12.5 


14.3 


15.3 


T52 


11 12 27.72 


-76 44 22.3 


6.84 


G9 


5410 


2.23 


31 






0.110 


8.2 


9.9 


11.4 


13.0 


T54 e 


11 12 42.69 


-77 22 23.1 


7.88 


G8 


5520 


2.41 


30 


30* 




0.083 


10.2 


12.4 


13.8 


15.3 


T33A+B e 


11 08 15.10 


-77 33 53.2 


6.18 


G7 


5630 


2.51 


30 


30 




0.086 


8.4 


10.4 


11.8 


12.8 


T21 


11 06 15.41 


-77 21 56.8 


6.42 


G5 


5770 


2.54 


30 


30 




0.108 


7.8 


9.7 


11.2 


12.1 


T26 e 


1 1 07 20.74 


-77 38 07.3 


6.22 


G2 


5860 


2.57 


90' 


45J 


2 


0.072 


8.8 


10.9 


12.6 


14.3 


HD 93828 


10 46 37.95 


-77 36 03.6 


7.48 


FO 


7200 


2.87 


31 


31 


2 


0.082 


9.9 


12.1 


13.7 


15.5 


T41 e 


1 1 09 50.03 


-76 36 47.7 


7.15 


B9 


10500 


> 3 


31 


31 


2 


0.094 


8.9 


11.0 


12.5 


13.6 


HD 96675 


11 05 57.81 


-76 07 48.9 


6.97 


B6.5 


13500 


> 3 


31 


31 


2 


0.099 


8.5 


10.4 


11.6 


12.6 


T32 


11 08 03.30 


-77 39 17.4 


5.94 


B9.5 


10010 


>3 


90' 


31 


2 


0.072 


8.4 


10.7 


12.5 


14.0 



a From Luhman (2004). 

b Unless stated otherwise, estimate based on T c n and the models of D'Antona & Mazzitelli (1997); see text for more detail. 
c Dichroic used: (1) N90C10 or (2) VIS. 
d In AT-band. 

e Object was known as a binary from previous studies. 

f Estimate based on absolute AT-band magnitude and the models of D'Antona & Mazzitelli (1997); see §3.2 for more detail. 
g Object was known as a triple from previous studies. 
h Complete name is 2MASS II 1 103481-7722053. 
1 This observation was done in the NB2.17 filter. 
J This observation was done in the NB1.64 filter. 

^ We did not detect the substellar companion to this source reported by Luhman et al. (2006) as it is below our detection limits. 

The spectroscopic data were reduced in a standard way, a parametric PSF allowed to vary with wavelength. Compar- 

except that an unfortunate ghost at the location of the com- isons with spectra extracted by simply summing up all flux in 

panion position in the data of CHXR 18N data prevented us a large aperture yielded consistent results. Wavelength cal- 

from using the B position as sky for A. Instead, we fit a lin- ibration was made using night sky OH emission lines and 

ear function to the sky along the slit. This procedure did not the atlas of Rousselot et al. (2000). The spectra were cor- 

leave any visible residuals from the sky OH lines. The spectra rected for extinction, using the R = 3.1 parametric formula 

were extracted using an optimal extraction algorithm, fitting from Fitzpatrick (1999) and assuming the same extinction as 
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measured for the primaries by Luhman (2004, Ay=0.60 for 
T54 and A j=0. 14 for CHXR 1 8N). 

The motivation for choosing the 172 mas wide slit was the 
need to ensure that the companion candidates would be lo- 
cated in the slit using blind acquisitions, given their relatively 
large contrasts and separations with respect to their primaries. 
At a separation of ~2.5", a small uncertainty of ~1° in po- 
sition angle would translate into a ^50 mas uncertainty in 
position, comparable to the PSF FWHM. In addition, as the 
AO correction varies strongly with wavelength, the resulting 
PSF is wavelength-dependent; this effect was amplified by the 
relatively high airmass (1.7) of our observations, which fol- 
lows from the southern declination of the targets (-77°). The 
wide slit used thus considerably simplifies the wavelength- 
dependent relative flux calibration, but also adds the compli- 
cation that the spectral resolution becomes PSF dependent, 
in our case R = 200-500, and with a complicated line-spread 
function. This limits the accuracy to which we can remove 
telluric absorption lines by dividing by the calibration stars 
(and multiplying by a model), as the spectral resolution for 
the science and calibration observations differ. To mitigate 
this effect, we convolved the calibration spectra with a Gaus- 
sian until the telluric residuals were minimized (although not 
entirely removed). Since our main interest lies in the shape of 
the spectrum, however, and not in individual absorption lines, 
we did not try to optimize further the telluric feature removal. 

3. ANALYSIS 
3.1. Identification and measurements of multiple systems 

Identification of multiple systems was first done by visual 
inspection of the images. Then, for well-separated multiple 
systems (> 0.5"), subtraction of the stellar signal using the 
PSF of another component was performed to search for addi- 
tional tight, low-mass companions to each of the components. 
While this approach led to very small residuals, no additional 
companion was found. For apparently single stars, PSF sub- 
traction using the images of other targets was attempted but 
yielded little improvement due to the large variability of the 
PSF throughout the three nights and again, no additional com- 
panion was found. Although the PSFs of a few targets are 
elongated, which could indicate unresolved binarity, it is diffi- 
cult to make concrete assertions as there are several instances 
of intrinsically elongated PSFs in our data (i.e. elongation 
visible in both components of a binary system). 

The angular separation, position angle, and magnitude dif- 
ference between the components of multiple systems were de- 
termined using a custom PSF extraction/subtraction IDL rou- 
tine. For well-separated systems, this routine determines the 
parameters that minimize the residual noise after subtraction 
of one component from another. For close-separation sys- 
tems, the routine first reconstructs the PSF, using the then best 
guess of the separation and flux ratio of the components, from 
the parts of each component's PSF that are uncontaminated 
by the other component, and then determines the parameters 
that minimize the residual noise after subtraction of the re- 
constructed PSF from the image. This process is iterated until 
convergence is reached. The random errors on each parameter 
were estimated from the dispersion of the measurements made 
on the six individual images of each sequence. By compar- 
ing the FWHM and encircled energy of the primary stars and 
their companions as a function of angular separation, we did 
not find evidence for significant anisoplanatism effects below 
6", and thus we did not apply any correction to the measured 



photometry for such effects. 

3.2. Mass estimates 

The masses of all targets were estimated from their effec- 
tive temperatures and the predictions of pre-main sequence 
evolution models; this approach was preferred over a deter- 
mination based on luminosity as it is less affected by unre- 
solved binarity or uncertainties in extinction or excess emis- 
sion. Choosing which pre-main-sequence evolution model to 
use is a delicate task as these models are highly uncertain at 
young ages and there are significant differences between the 
different ones available. For the present study, as one of its 
goals is to investigate the variations of the properties of mul- 
tiple systems as a function of mass, the masses should ideally 
be estimated in a consistent manner across the entire range 
of masses observed (^0.02^1 M Q for primaries and compan- 
ions). We have thus opted for the models of D'Antona & 
Mazzitelli (1997) 4 , whose wide mass coverage (0.017-3 Mq) 
enables determination of the masses of all but three of the pri- 
mary stars observed and all but one of the companions found 
(see §4). For completeness, a discussion of the effects of us- 
ing other evolution models to determine masses is presented 
in §A. For three targets for which the spectral type (hence ef- 
fective temperature) is unknown, the masses were determined 
based on absolute TT-band magnitude. To do so, the lumi- 
nosities from the models of D'Antona & Mazzitelli (1997) 
were converted to K magnitudes using the bolometric correc- 
tions and colors tabulated in Kenyon & Hartmann (1995) for 
TeB > 3000, and in Bessell (1991) for r eff < 3000. The esti- 
mated masses of all targets are indicated in Table 1 . An age 
of 2 Myr (Luhman 2004) was assumed for all targets to pro- 
duce these estimates. This simple assumption introduces ad- 
ditional uncertainties as in reality there is likely an age spread 
among the Cha I members (~l-5 Myr, see Fig. 15 of Luh- 
man 2004). By comparing the masses obtained for ages of 1 
or 5 Myr to those obtained for an age of 2 Myr using the mod- 
els of D'Antona & Mazzitelli (1997), we estimate that these 
uncertainties are of the order of 10-15%. 

The mass ratios of the secondary (and tertiary) components 
of multiple systems were estimated based on their contrast 
with respect to their primary. More precisely, the observed 
TT-band contrast of the companion was added to the model K- 
band magnitude of the primary (i.e. based on its effective tem- 
perature), and the corresponding companion mass was then 
obtained from the models; this approach is less sensitive to 
uncertainties in distance and extinction than one relying on 
observed absolute magnitudes. However, the observed flux 
ratios may be slightly different from the true photospheric flux 
ratios if there is excess emission from circumstellar material 
around one or both components of a system. For reference, 
we point out that an error of 0.2 mag on the photospheric flux 
ratio would typically lead to a difference of roughly 15% in 
mass ratio. As discussed before for the primary mass, the 
uncertainties on the ages of the systems also introduce uncer- 
tainties on the estimated mass ratios. We estimate that these 
uncertainties are of the order of 15-20%, using the evolu- 
tion models of D'Antona & Mazzitelli (1997). For the ter- 
tiary component of CHXR 68 and the secondary component 
of Hn 21 W, we have used the effective temperatures (3524 K 
and 3024 K, resp.) reported in Luhman (2004) to derive their 
masses. 

4 We have used the 1998 update of these models for masses in the range 
0.017-0.9 M Q . 
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3.3. Detection limits 

The detection limits as a function of angular separation 
from the target stars were obtained by calculating the resid- 
ual noise in annuli centered on the target after subtraction of 
an azimuthally symmetric PSF profile and convolution of the 
residual image by a circular aperture of radius 0.04". The de- 
tection limits were determined as 5 times the residual noise. 
For multiple systems, the secondary and/or tertiary compo- 
nent was subtracted from the image using an analytic PSF 
model prior to the calculation of the detection limits. Since 
at the smallest separations the variance of the residual noise 
may be underestimated by this procedure as there are too few 
independent spatial elements within the corresponding annuli, 
we have required that the detection limits within a radius of 2 
FWHM be no better than the PSF intensity profile. As stated 
before, we found no evidence for significant anisoplanatism 
within the angular separation range probed, and thus we did 
not apply any correction to the detection limits for such an 
effect. 

The /T-band detection limits are indicated in Table 1 for four 
angular separations. We provide detection limits at 0.1" even 
for targets for which the PSF FWHM exceeds 0.1" since for 
these targets the PSF is very smooth and very finely sampled 
by the 0.013" pixels such that companions with separations 
slightly below one FWHM could be detected, provided they 
are above the primary PSF intensity profile. The median K- 
band detection limits are 12.2 mag at 0.1", 14.1 mag at 0.25", 
14.8 mag at 0.5", and 15.0 mag at 1" and beyond. In terms 
of contrast, the median limits are 2.2, 4.4, 5.1, and 5.3 mag, 
respectively. Mass ratio limits have been calculated for each 
target based on their K s contrast limit in the manner described 
above. For each of four mass bins and two angular separa- 
tions, the minimum mass ratios that could have been detected 
for 50% (^50%) and 90% (<?9o%) of the targets are reported in 
Table 2; the latter value can be used as an estimate of our 
completeness limit. The higher mass ratio limits at small sep- 
arations for the last mass bin of Table 2 arise from a combi- 
nation of higher speckle noise for these brighter stars and a 
brightening of ^2-3 M Q stars at ^2 Myr as they settle on the 
zero-age main sequence and the CNO cycle reaches equilib- 
rium (L. Siess 2007, private communication). 



TABLE 2 
Mass ratio limits 


Mass bin 


Mass ratio limit (qso%/l90%) 


(Mq) 


0.1" 0.25" 


0.5" 


1" 


0.10-0.23 
0.23-0.55 
0.55-1.28 
1.28-3.00 


0.23/0.44 0.12/0.24 
0.13/0.20 0.04/0.08 
0.11/0.23 0.04/0.06 
0.29/0.60 0.04/0.14 


0.09/0.21 
0.03/0.05 
0.02/0.04 
0.02/0.03 


0.09/0.21 
0.03/0.06 
0.02/0.03 
0.01/0.02 



4. RESULTS 

4. 1 . Resolved multiple systems 

We have identified a total of 30 candidate binary systems 
and 7 candidate triple systems among the 126 targets of our 
sample. No higher-order multiple system was found. The ma- 
jority of the binaries and all but two (T31 and T39) of the 
triples are new discoveries. The triple systems CHXR 28, 
CHXR 68A and T26 were previously known binaries for 
which one of the components was itself resolved into a tight 
binary. The properties of the binary and triple systems ob- 
served are summarized in Tables 3 and 4. Only the random 
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0.1 1.0 10.0 

Angular separation (") 

FIG. 1. — Separation and contrast, relative to their primary star, of all the 
companions identified in this study. The two curves indicate detection limits 
that were obtained for at least 90% {top) and 50% (bottom) of the targets. 

measurement errors are reported in Table 3, a systematic un- 
certainty of 0.226% in separation and 0.5° in position angle 
should be added (Masciadri et al. 2003). The separation and 
K s contrast of all companions identified are shown in Figure 1 . 
Figure 2 presents images of the tightest binary found in our 
survey, 2M1 110-7722, which is marginally resolved with a 
separation of 0.059" and contrast of 0.69 mag in K s . Images 
of the triples systems identified are shown in Figure 3. The 
faint companion candidate around the binary T54 and the faint 
companion candidate around CHXR 18N are omitted from ta- 
ble 4 and figures 1, 3 as follow-up observations have indicated 
that they are background stars (see §4.3). They are however 
listed in table 3 for completeness. 

The triple systems were named based on hierarchical lev- 
els. The components of the wide pairs were assigned the cap- 
ital letters A,B and the components of the tight pairs were as- 
signed the lower case letters a,b; the brighter component of 
each pair bears the letter A or a. The configuration and prop- 
erties of the triple systems, split into tight and wide pairs, are 
indicated in Table 4. For the wide pair, the tight subsystem is 
considered as a single entity whose mass is equal to the sum of 
the masses of its components, and the separation of the wide 
system is calculated with respect to the center of mass of the 
tight subsystem. Among our sample of triple systems, the two 



• 






• 




■ % 



FIG. 2. — K s - (top) and //-band (bottom) AO images of the tight bi- 
nary system 2M1 110-7722 (first column). The other columns show the 
same images after subtraction of the best-fit single-component (second col- 
umn) and two-component (third column) analytic PSF models consisting of 
Gaussian+Moffat functions; the display intensity scale is ±10% of the pri- 
mary PSF peak. The PSF models are allowed to be asymmetric. Each panel is 
0.5" on a side. The smaller residuals after subtraction of the two-component 
PSF model clearly indicate that this is a binary system rather than a single 
star with an elongated PSF. 
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TABLE 3 

Measurements and properties of multiple systems 



Name 


Separation" 


Position Angle b 


AK S 


AH 




P' 




(") 


(°) 


(mag) 


(mag) 




(yr) 


Binaries 


2M1 110-7722 


0.059 ±0.002 


30.1 ±1.5 


0.69 ± 0.04 


0.57 ±0.19 


0.46 


50 


B53 


0.295 ±0.001 


235.2 ±0.1 


1.52 ±0.03 


1.55 ±0.01 


0.24 


600 


CHXR 15 


0.304 ±0.001 


199.9 ±0.2 


0.65 ±0.01 


0.63 ±0.01 


0.53 


700 


CHXR 26 


1.396 ±0.001 


33.8±0.1 


0.47 ±0.01 


0.45 ±0.01 


0.51 


6000 


CHXR 30A 


0.459 ±0.001 


23.2±0.1 


0.85 ±0.01 


0.95 ±0.02 


0.41 


700 


CHXR 37 


0.079 ±0.001 


81.0±1.1 


1.07 ±0.01 


1.09 ±0.04 


0.32 


50 


CHXR 40 


0.151 ±0.002 


66.1 ±0.6 


0.12±0.03 


0.19±0.03 


0.89 


130 


CHXR 47 


0.175 ±0.001 


334.8 ±0.2 


0.47 ±0.05 


0.46 ±0.01 


0.73 


120 


CHXR 49NE 


0.267 ±0.001 


53.2±0.1 


0.05 ±0.01 


0.04 ±0.01 


0.96 


300 


CHXR 59 


0.148 ±0.001 


347.4 ±0.2 


0.02 ±0.01 


0.06 ±0.03 


0.97 


160 


CHXR 62 


0.120 ±0.002 


135.1 ±1.0 


0.06 ±0.09 


0.12 ±0.09 


0.93 


130 


CHXR 71 


0.572 ±0.001 


70.0 ±0.1 


1.63 ±0.01 


1.75±0.01 


0.23 


1600 


CHXR 79 


0.885 ±0.001 


211.2±0.1 


2.45 ± 0.02 


2.10±0.03 


0.11 


2000 


Cha Ha 2 


0.167 ±0.001 


40.4 ±0.2 


0.11 ±0.02 




0.77 


300 


Hn 13 


0.132±0.001 


318.2±0.2 


0.09 ±0.01 


0.17±0.01 


0.81 


190 


Hn21W 


5.495 ±0.004 


69.2 ±0.1 


0.95 ±0.02 


0.91 ±0.02 


0.71 


40000 


Hn4 


0.211 ±0.001 


296.1 ±0.3 


0.04 ±0.01 


0.02 ±0.18 


0.95 


300 


ISO 126 


0.292 ±0.001 


232.4 ±0.2 


0.66 ±0.01 


0.49 ±0.01 


0.46 


400 


T21 


0.140 ±0.006 


126.1 ±2.3 


2.16±0.12 


2.05 ±0.11 


0.33 


60 


T27 


0.787 ±0.001 


13.5 ±0.1 


0.67 ±0.01 


0.50 ±0.01 


0.46 


1800 


T'J 
1 J 


7 7 i « -t- n nn7 
Z.ZIO It u.uuz 


7so nidi 
zoy.u it u. l 


l.oU It U.U1 


i 1 7 4- n m 
1 . t / It u.uz 


u. 10 


vuuu 


T33A+B 


2.434 ±0.001 


284.8 ±0.1 


1.93 ±0.04 


0.66±0.12 


0.38 


4000 


T41 


0.788 ±0.001 


326.7 ±0.1 


3.05 ±0.03 


3.33 ±0.06 






T43 


0.796 ±0.001 


352.2±0.1 


1.41 ±0.01 


1.43 ±0.01 


0.22 


2000 


T45 


0.752 ±0.006 


52.2 ±0.4 


2.67 ±0.02 


2.34 ±0.03 


0.09 


1700 


T46 


0.123 ±0.001 


241. 6 ±0.6 


1.57 ±0.02 


1.38 ±0.08 


0.18 


110 


T5 


0.159±0.001 


342.9 ±0.1 


0.05 ±0.01 


0.04 ±0.01 


0.95 


190 


1 

1 Jl 


1 Q77 4_ n nm 


i ^7 ^ 4- n 1 

lOZ. J It U. 1 


7 ^ 4- n n i 
Z.7J It u.ui 


7 34 4- n ni 

Z. 74 It U.U7 


U. 1U 


uUUU 


T54 d 


n 747 _i_ n nm 


746 «J-m1 
iHU. J zn u.z 


l 45 4- n 04 
i .tj it u.w 


1 54 4- n (11 

l.Jt It \J.\J 1 


0.52 


130 


T6 


5.122±0.003 


122.9 ±0.1 


3.79 ±0.03 


3.90 ±0.06 


0.04 


16000 






Triples 
























V n A K Zo /\a,/\D 


c\ ■\A r i nni 

U. 14J ± U.UUl 


356.6 ±1.0 


0.40 ±0.05 


U.4U ± U.U4 


n 73 

0. 1 J 


1 fin 
tuu 


CHXR 28 Aa,B 


1.818±0.003 


115.9 ±0.1 


0.32 ±0.04 


0.31 ±0.10 


0.77 




CHXR 68A Aa,Ab 


0.101 ±0.001 


13.6±1.5 


0.22 ±0.02 


0.33 ±0.07 


0.84 


60 


CHXR 68A Aa,B 


4.367 ± 0.005 


213.0±0.1 


0.95 ±0.12 


0.27 ±0.57 


0.58 




CHXR 9C A,Ba 


0.852 ±0.001 


81.3 ±0.1 


0.36 ±0.04 


0.45 ±0.11 


0.63 




CHXR 9C Ba,Bb 


0.1 30 ±0.002 


76.1 ±0.4 


0.70 ±0.04 


0.67 ±0.11 


0.46 


170 


T26 A,Ba 


4.557 ±0.007 


202.1 ±0.1 


5.00 ±0.03 


4.55±0.12 


0.03 




T26 Ba,Bb 


0.066 ±0.005 


177.9 ±1.8 


0.06 ±0.32 


0.17±0.16 


0.97 


90 


T31 A,Ba 


0.652 ±0.001 


178.7 ±0.1 


1.74 ±0.01 


1.22 ±0.01 


0.16 




T31 Ba,Bb 


0.100±0.001 


231.0±0.3 


0.61 ±0.02 


0.65 ±0.07 


0.7 


170 


T39 Aa,Ab 


1.242 ±0.001 


19.2±0.1 


0.31 ±0.01 


0.31 ±0.01 


0.67 


4000 


T39 Aa,B 


4.497 ±0.001 


70.8 ±0.1 


0.77 ±0.02 


0.87 ±0.02 


0.34 








.Unrelated hackpround stars 






















T54 A,c d 


2.383 ±0.007 


81.10±0.15 


7.06 ±0.15 


7.13±0.15 






CHX 18N A,b d 


2.738 ±0.006 


356.3 ±0.1 


5.40 ±0.03 


5.25 ±0.03 







a Random uncertainties are listed. An additional systematic uncertainty of 0.226% is present. 

b Random uncertainties are listed. An additional systematic uncertainty of 0.5° is present. 

c Estimate based on AK S using the models of D'Antona & Mazzitelli (1997), see §3.2 for more detail. 

d AH and A^ s shown were measured in 2007 follow-up imaging data. Contrasts of AJ = 1 .65 ± 0.06, AJ = 7.6 ± 0.3, and AJ = 5.28 ± 0.07 were further 
measured for T54 B, T54 c, and CHXR 18N b, respectively. 

TABLE 4 
Properties of triple systems 



. Tight subsystem . . . . Wide subsystem . 



Name 


Configuration 


OA") 




PI (yr) 








kip; 


CHXR 28 


(Aa,Ab),B 


0.143 


0.73 


100 


1.85 


0.45 


4000 


40 


CHXR 68A 


(Aa,Ab),B 


0.101 


0.84 


60 


4.41 


0.32 


16000 


260 


CHXR 9C 


A,(Ba,Bb) 


0.130 


0.46 


170 


0.89 


0.92 


2000 


12 


T26 


A,(Ba,Bb) 


0.066 


0.97 


90 


4.59 


0.05 


12000 


130 


T31 


A,(Ba,Bb) 


0.100 


0.70 


170 


0.68 


0.27 


1400 


8.2 


T39 


(Aa,Ab),B 


1.242 


0.67 


4000 


4.21 


0.20 


21000 


5.3 
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FIG. 3. — Triples systems identified in our survey. The intensity scale is linear, except for T26 which is shown on a logarithmic intensity scale. 



possible configurations, (Aa,Ab),B and A,(Ba,Bt>), are ob- 
served three times each. 

A pseudo orbital period P' was calculated for each system 
assuming that the actual orbital semi-major axes are equal to 
the projected separations observed. The mass ratios of the 
companions and the pseudo orbital periods are indicated in 
Tables 3 and 4. Three binary systems have a pseudo-orbital 
period of ^50 yr and thus for these systems, depending on 
their orbital phase and orientation, a partial orbit could be 
measured within the next 10-20 years to derive dynamical 
constraints on their masses. The tight pair of CHXR 68A has 



the shortest pseudo-period among the triple systems observed 
and could allow a partial orbit determination within 10-20 
years. 

4.2. Probability of chance alignment 

We have used star counts based on the Two Micron All-Sky 
Survey (2MASS) Point Source Catalog (PSC) to estimate the 
probabilities that the candidate companions identified in our 
images are unrelated background stars rather than real com- 
panions; a similar approach was used by Correia et al. (2006). 
For each target in our sample, we have retrieved all sources 
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from the 2MASS PSC that lie within a radius of 15' and we 
have built a cumulative distribution of the number of sources 
as a function of K s magnitude. These cumulative distributions 
were then converted to surface densities. For most of our tar- 
gets these surface densities are slight overestimates as their 
immediate vicinity is usually less densely populated, owing to 
higher extinction, than their broader surroundings extending 
to 15'. Based on these surface densities, we have calculated, 
for each of the 44 candidate companions identified, the prob- 
ability that a background star at least as bright as the compan- 
ion and within the same angular separation would have been 
found, and detected, around each of the 126 targets. For the 
vast majority of candidates and targets, these probabilities are 
below 10~ 4 . Then for each candidate companion, by summing 
these background source detection probabilities over the 126 
target stars, we obtained the expected total number of back- 
ground sources of the same brightness or brighter and within 
the same angular separation that should have been detected by 
our survey. 

For all but 5 of the candidate companions, this expected 
number is at the level of 0.03 or lower, indicating that the 
candidate companions are almost certainly bound to their pri- 
mary. The 5 candidate companions with higher expected 
numbers are the faint tertiary around T54 (expected number 
of equivalent background sources detected of 0.5), the faint 
source around CHXR 18N (0.4), the source around T6 (0.25), 
the tight binary system around T26 (0.25), and the most dis- 
tant source around T39 (0.08). Given these expected numbers, 
denoted A, and assuming Poisson statistics, the probabilities 
that the distant companions are truly bound objects are sim- 
ply given by e~ x , i.e. the probability to detect zero equiv- 
alent background source. For the above candidate compan- 
ions, these probabilities are 0.60, 0.67, 0.78, 0.78, and 0.92, 
respectively. The five sources with the lowest probabilities 
of being bound are discussed in more detail in the next sec- 
tion; all other candidates are considered bound for the present 
study. 

The above analysis for the binary systems 
CHXR 68A/CHXR 68B and Hn 21W/Hn 21E gives ex- 
pected numbers of equivalent background stars of 0.08 and 
0.4, respectively. However, as both components of these 
systems are confirmed members of Cha I and the surface 



density of members of Cha I is much smaller than the surface 
density of foreground/background stars, the above analysis 
does not apply to these systems and we consider them as 
almost certainly physically bound. 



4.3. 



Follow-up of the widest companion candidates 
section additional astrometric and snectrosc 



In this section, additional astrometric and spectroscopic 
data are used to establish the nature of the five candidate com- 
panions with the lowest probabilities of being bound based on 
stellar surface density statistics. This analysis is based on the 
imaging and spectroscopic follow-up of the sources T54 and 
CHXR 18N that was mentioned previously, as well as on data 
retrieved from the HST archive and from the literature. 

We have retrieved and analyzed HSTAYFPC2 images for 
the stars CHXR 18N, T6 and T26. In some of those images, 
the peak of the primary star PSF is saturated and thus cannot 
be used for centroid determination. In those cases, the cen- 
troid of the primary star was found by maximizing the cross- 
correlation of the secondary mirror support diffraction spikes 
after a 180° rotation of the image. The candidate compan- 
ion around CHXR 18N is also saturated in the HST/WFPC2 
images, but its diffraction spikes are too faint to be used for 
PSF registration; in this case, the centroid was determined by 
cross-correlation, within an annulus excluding the saturated 
pixels, with the PSF of an unsaturated nearby star in the im- 
age. In all other cases, i.e. for unsaturated PSFs, the centroid 
was determined by fitting a 2-D Gaussian function. For this 
analysis, the tight binary around T26 was treated as a sin- 
gle object and measurements were made relative to its center 
of light. While not resolved in the HST images, the PSF of 
this tight binary shows hints of elongation. For the candidate 
companion T39 B, we have used relative astrometry measure- 
ments reported in Correia et al. (2006), based on VLT/NACO 
observations. All of the astrometric measurements are sum- 
marized in Table 5. 

Figure 4 shows a comparison of the measured separations 
and position angles of the five candidate companions with 
those expected for stationary background stars based on the 
proper motions of their primary. The proper motions of the 
primary stars were taken from the UCAC2 catalog (Zacharias 
et al. 2004), except for the star T54 for which it was taken 



TABLE 5 

Multiple epoch astrometry of the widest candidate companions 



System 


Epoch 


Separation" 


Position Angle" 


Source 




(yr) 


(") 


(°) 




T54 A,c 


2006.23098 


2.383 ± 0.007 ± 0.005 


81.10±0.15±0.5 


1 




2007.36340 


2.409 ± 0.003 ± 0.005 


80.62 ± 0.07 ±0.2 b 


1 


CHXR 18N A,b 


2000.21481 


2.76 ±0.03 


356.3 ±0.7 


2 




2006.23382 


2.738 ± 0.006 ± 0.006 


356.3±0.1±0.5 


1 




2007.46197 


2.747 ± 0.003 ± 0.006 


356.12±0.06±0.5 


1 


T6 A,B 


1998.78593 


5.118±0.012 


123.00±0.15 


2 




2006.23380 


5.122 ±0.003 ±0.011 


122.9 ±0.1 ±0.5 


1 


T26 A,(Ba,Bb) c 


1997.92682 


4.63 ±0.03 


202.3 ±0.4 


2 




2000.09098 


4.589 ±0.012 


202.30 ±0.15 


2 




2006.23354 


4.587 ±0.010 ±0.010 


202.0 ±0.1 ±0.5 


1 


T39 Aa,B 


2003.05233 


4.504 ±0.002 ±0.010 


70.8 ±0.1 ±0.5 


3 




2006.22815 


4.497 ±0.001 ±0.010 


70.8 ±0.1 ±0.5 


1 



REFERENCES. — (1) This work; (2) HSTAVFPC2 archive; (3) Correia et al. 2006. 

a When two uncertainties are quoted, the first is the measurement uncertainty and the second is the systematic uncertainty. 

b The field of view orientation of the 2007 images was calibrated relative to the 2006 images using the inner binary of this system; a correction of -0.6° was 
applied to the second epoch orientation obtained from the image header. The systematic uncertainty on the position angle of the second epoch measurements 
thus corresponds to the random uncertainty of the position angle of the inner binary. 

c Relative position for the center of light of component Ba,Bb. 
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FIG. 4.— Relative positions of the candidate companions T54 c (a), CHXR 18N b (b), T6 B (c), T26 BaBb (rf) and T39 B (e) as a function of time. Open 
diamonds mark the observed separations (top) and position angles {bottom) of the candidate companions at the epochs for which we have measurements. The 
solid lines indicate the expected separation and position angle of a stationary background source as a function of time; the dotted lines reflect the uncertainties on 
the proper motion and distance of the primary stars. 
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FIG. 5. — //-band spectra of the companion candidates CHXR 18N b 
and T54 c. The spectra are normalized at 1.6 fim and offset by 0.15 for 
clarity. The middle spectrum is an observed Ml. 5 V comparison spectrum 
(HD 36395) from the IRTF spectral library (Cushing et al. 2005), smoofhened 
to the approximate resolution of the other spectra. 

from Frink et al. (1998). 5 Figure 4 clearly indicates that 
T6 B, T26 BaBb and T39 B are truly bound to their pri- 
maries, and we consider them as such for the rest of the anal- 
ysis. However for the candidates T54 c and CHXR 18N b, 
the comparison is rather inconclusive. For T54 c, the re- 
duced x values for the bound and not bound scenarios are 
6.7 and 4.5, respectively, making it impossible to differentiate 
between them. For CHXR 18N, the reduced \ 2 values are 1.2 



and 6.4, respectively, revealing a better agreement with true 
physical association but with a low level of significance. 

The spectra of the companion candidates around T54 and 
CHXR 18N are shown in Figure 5. If the companion candi- 
dates had been real companions, their K-band contrasts would 
have implied late-M spectral types and masses of ~0.025 M Q 
and ^0.015 M@, respectively. Late M-type objects have a 
characteristic triangular-shaped //-band spectrum (Cushing 
et al. 2005) that is further emphasized by the lower surface 
gravity of young objects (Lodieu et al. 2008). We com- 
pared the spectra of the two companion candidates to the 
DUSTY models by Chabrier et al. (2000) and to observed 
near-infrared spectra from the IRTF spectral library (Cush- 
ing et al. 2005), and found the spectra of these two compan- 
ion candidates to be clearly incompatible with spectral types 
later than M5. The best-fit spectrum from the IRTF library 
is also shown in Figure 5, and comes from the Ml. 5 V star 
HD 36395. We conclude that the companion candidates T54 c 
and CHXR 18N b are likely background dwarfs, located at 
^350 pc and ^800 pc, respectively. 

5. DISCUSSION 

5.1. Overall multiplicity frequency 

The multiplicity fraction (MF) and the mean number of 
companions per primary star, or companion star fraction 
(CSF), are indicated in Table 6 for our raw and complete 
Cha I samples. For our complete sample, defined as includ- 
ing only targets for which we have reached a contrast limit at 
least as good as the top curve (90%) of figure 1 and counting 
only companions falling above this curve and in the separation 
range 0.1"-6", the MF is 0.27!°$ and the CSF is 0.32+ 06 



--0.05- 



5 A significantly different value is reported for T54 in the Tycho-2 catalog 
(H0g et al. 2000), but we did not use this value as Skiff (2007) pointed out in 



a note that the astrometry of this source in the Tycho-2 catalog was wrong. 
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TABLE 6 

Multiplicity fractions of various samples 



Study 



Region S B T MF 



CSF 



h 



This work 



29+ 05 



33+ - 06 



17 +008 



27+ 05 
u - -0.04 



a -19+0.06 
u - ,z -0.05 



17+0O9 
u - -0.07 



Raw sample . . 

Cha I 90 30 6 

Complete sample" 

This work Cha I 80 25 5 

Restricted sample: 0.13" < 6 < 6.4", AK < 2.5, and K < 10.4 

This work Cha I 48 19 2 0.30*>°<j 0.33^ ^ 

Ratzka et al. (2005) p Oph 103 34 2 0.26+™ 0.27;™$ 

Kohler et al. (2000) Sco-Cen 59 27 2 0.33+" ; °| 0.35;*<" 

Kohler & Leinert (1998) Taurus 50 19 1 0.29™° °- 30 3)06 
Combined sample ■■■ 260 99 7 0.29^ ; || OJl^jf 

Restricted sample: 0.1 < M/M Q < 2.0, 60 < 0/AU < 500, and AK < 3 

This work Cha I 83 13 0.14^" 

Kohler et al. (2006) ONC 105 5.4 0.05™™ 

Restricted sample: 32 < fl/AU < 640 and q > 0.1 

This work Cha I 58 18 0.24^ ^ 

Duchene et al. (1999) IC 348 66 8.5 0.13 



10 +0 " 
06+o o" 

07+ 08 
()O5 +0.f3 

-QQ4 
07+0 03 
u - -0.02 
5 



-0 04 



a Considering only the 1 10 targets for which we have reached a contrast limit at least as good as the top curve (90%) of figure 1 and counting only companions 
falling above this curve and in the separation range 0.1 "-6". 



We first compare the MF and CSF in Cha I and in three 
other dispersed star-forming regions using the results of ear- 
lier surveys. However, as previous surveys have typically tar- 
geted brighter stars, had lower angular resolution, and were 
less sensitive, a direct comparison is impossible and we are 
forced to consider a more limited subset of our observations. 
Ratzka et al. (2005) conducted a multiplicity survey of 158 
members of the p Ophiuchus dark cloud using speckle imag- 
ing on a 3.5m telescope; all of their targets have a K-bwd 
magnitude <10.4. Here we consider only their restricted sam- 
ple, defined as comprising only those targets with high likeli- 
hood of membership in p Oph and with flux ratio limits <0.1 
at >0.13", and only those companions in the separation range 
0.13"-6.4" and TT-band flux ratio >0.1. Kohler et al. (2000) 
conducted a multiplicity survey of 118 X-ray selected stars 
in the Scorpius-Centaurus OB association also using speckle 
imaging on a 3.5m telescope. Again we consider only their 
restricted sample, defined as above and with a similar range 
of target magnitudes. Finally, Kohler & Leinert (1998) ob- 
tained similar observations of the Taurus star-forming region. 
Although they do not provide an equivalent restricted sample 
in their paper, we have constructed one using their Tables 1 
and 2. The statistics of these three samples are shown in Ta- 
ble 6, along with the results for our Cha I sample after apply- 
ing the same restrictions. As it can be seen in this table, the 
MF values are very similar in all four regions: the probabili- 
ties that the MF in Cha I and in p Oph, Sco-Cen, and Taurus 
were drawn from the same underlying binomial distribution 
are 0.51, 0.86, and 0.85, respectively. 6 Thus it appears that 
for brighter targets (typically earlier than ^M2-M3) and in 
the separation range 20-1000 AU, the multiplicity fractions 
are the same in most dispersed star-forming regions. Given 
that the four samples considered above are statistically indis- 
tinguishable, it is legitimate to combine them to obtain a more 
statistically significant sample of 257 singles, 97 doubles, and 
8 triples; this yields MF=0.290^;^ and CSF=0.31^. 

Previous multiplicity studies have revealed that the MF in 
dense star-forming clusters is significantly below that in more 
dispersed regions; we investigate this here by comparing our 
results with observations of the Orion Nebula Cluster (ONC) 

6 Calculated using the "exact" hypothesis test for comparison of two bino- 
mial distributions detailed in Appendix B2 of Brandeker et al. (2006). 



and IC 348. Kohler et al. (2006) surveyed 228 stars in the 
periphery of the ONC using adaptive optics, and they com- 
plemented their study with the speckle imaging results of Petr 
et al. (1998) for the cluster core. Given that these authors 
found only a small and statistically not very significant differ- 
ence in binarity between the cluster core and its periphery, we 
consider only their combined sample that incorporates stars in 
both the core and the periphery. For stars with high likelihood 
of membership in the cluster and with masses in 0.1-2 M Q , 
for the separation range 60-500 AU, and for contrast below 
roughly 3.5 mag in K, they report an MF of 0.051 ±0.027. 
Applying these restrictions to our Cha I sample, we obtain 
MF=0.14^. The probability that both MFs were drawn 
from the same distribution is only 0.04; the overabundance 
of binaries in Cha I relative to the ONC, by a factor ~2.7, is 
thus statistically significant. Duchene et al. (1999) observed 
66 stars in the IC 348 cluster using adaptive optics; their target 
sample comprises stars of roughly the same masses as ours. 
For the orbital separation interval 32-640 AU and for mass 
ratios >0.1, they report 8.5 binaries after correction for in- 
completeness; this corresponds to MF=0. 13^04. For equiva- 
lent restrictions in our Cha I sample, we obtain MF=0.24±g {J 6 . 
Thus the MF in Cha I is also higher than that in IC 348, by a 
factor ^1.8; the probability that both MFs were drawn from 
the same distribution is 0.07. Our observations thus support 
earlier findings that the MF is lower in dense star-forming 
clusters than in dispersed star-forming regions. The binary 
statistics of the ONC, IC 348, and Cha I samples are summa- 
rized in Table 6. 

The ratio of the number of systems with at least three com- 
ponents to the number of systems with at least two compo- 
nents, denoted fa, is another interesting quantity to compare. 
This value is indicated in Table 6 for our Cha I sample as 
well as for the p Oph, Sco-Cen, and Taurus samples described 
above. The probabilities that the values of fa in Cha I and in 
p Oph, Sco-Cen and Taurus were drawn from the same distri- 
bution are all above 0.6; the differences are thus not statisti- 
cally significant. Combining the four samples, as done above 
for the MF, leads to a more statistically significant value of 
/ 3 = 0.08^. Correia et al. (2006) obtained AO imaging of 
58 known wide binaries (> 0.5") in a few star-forming re- 
gions. For their combined sample of 31 targets in Tau-Aur, 
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Cha I, and p Oph, they resolved 9 high-order multiples, lead- 
ing to fi w = 0.29^0 Qg. When considering only wide binaries 
in our sample (19 wide binaries, of which 6 have a third com- 
ponent), the corresponding value is / 3w = 0.32±g;['|, indistin- 
guishable from the value found by Correia et al. (2006). 

Comparison with the field population is slightly compli- 
cated as field multiplicity studies have typically concentrated 
their efforts on a narrow spectral type interval (e.g. Duquen- 
noy & Mayor (1991) for G stars; Delfosse et al. (2004) for M 
stars) and have explored a much wider orbital separation inter- 
val than our present observations. A meaningful comparison 
with the properties of our entire sample requires first taking 
the initial mass function (IMF) into account. This was done 
by Lada (2006) who reports an IMF-weighted MF of -0.33 
for field stars; the same value was also obtained by Reid & 
Gizis (1997) for the 8 pc sample. This number is based on 
statistics that should be complete over all orbital separations 
and for mass ratios >0.1. The restricted orbital separation 
range of our observations must then be taken into account. If 
we assume that the orbital separation distribution of all field 
multiple systems is the same as that reported for G-type stars 
by Duquennoy & Mayor (1991), then —41% of field multiple 
systems would have a separation within 16-1000 AU, and the 
corresponding field MF in that interval would be —0.14. Al- 
though this is a very rough comparison, it seems to indicate 
that within the 16-1000 AU interval, the MF in Cha I (0.27) 
is higher than in the field (—0.14), by a factor —1.85. The 
probability that the MFs in Cha I and in the field were drawn 
from the same distribution is below 0.04; in fact, it is likely to 
be even lower since the above assumption overestimates the 
number of M dwarf binaries in the 16-1000 AU range where 
studies reveal a deficit of wide companions among M dwarfs 
compared to G-type stars. 

5.2. Multiplicity as a function of mass 

To investigate the variation of multiplicity as a function of 
primary mass, we have divided the sample into 4 nearly-equal 
logarithmic mass bins and calculated the MF and CSF in each 
bin; the results are shown in Figure 6 and summarized in Ta- 
ble 7. There is a clear increase in multiplicity fraction with 
increasing primary mass: it rises from 0.16^ 05 to 0.63+Q23 
over the range —0. 1-3 M Q . This trend has been observed pre- 
viously both for other star-forming regions and for the field. 

We first compare these multiplicity fractions with values 
for field stars from the M-type stars survey of Delfosse et al. 
(2004); the visual multiplicity survey of 58 M2-M4.5 stars 
done by Fischer & Marcy (1992), counting only the compan- 
ions they found in the 10-1000 AU range; the multiplicity sur- 
vey of 41 M2-M4.5 stars done by Reid & Gizis (1997) using 
HST, sensitive to companions in the range 16-1000 AU; and 
the G-type stars survey of Duquennoy & Mayor (1991). As 
the surveys of Delfosse et al. (2004) and Duquennoy & Mayor 
(1991) are complete over all orbital separations, the MF val- 
ues they have reported must be corrected for the range of or- 
bital separations probed by our observations (—12-1000 AU) 
to make a meaningful comparison. Based on the orbital pe- 
riod distribution of Duquennoy & Mayor (1991), —44% of 
G-type multiple systems have a semi-major axis in the 12- 
1000 AU interval, while the same orbital separation range 
includes roughly 1/3 of the M dwarf multiple systems in 
Delfosse et al. (2004) (see their Fig. 4). We have thus mul- 
tiplied the MF values reported in Delfosse et al. (2004) and 
Duquennoy & Mayor (1991) by 0.33 and 0.44, respectively. 
The field values from these four studies are summarized in 
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FIG. 6. — Multiplicity fraction (solid histogram) and companion star frac- 
tion (dashed histogram) of our sample as function of primary mass. Open di- 
amonds show the MF of the numerical simulations of Goodwin et al. (2004a), 
adapted to the resolution limit of our observations (see §5.7 for detail). The 
cross symbols show field MF measurements corrected for the limited or- 
bital separation range probed by our survey (see §5.2 for detail); the values 
are taken from Delfosse et al. (2004) for the leftmost point, from Fischer 
& Marcy (1992) and Reid & Gizis (1997) for the middle point, and from 
Duquennoy & Mayor (1991) for the rightmost point. 



Table 7 and shown in Figure 6, which indicates that the in- 
crease of the MF as a function of primary mass in Cha I is 
similar to that of the field, although the values themselves are 
systematically higher in Cha I. 

A few recent high angular resolution multiplicity surveys of 
star-forming regions have targeted only low-mass stars, and 
their results are suitable for comparison with our results here. 
The Keck speckle interferometry survey of Konopacky et al. 
(2007) and the HST survey of Kraus et al. (2006) have targeted 
a combined sample of 15 Taurus members in the 0.1-0.2 Mq 
range and found 3 binaries within the separation range probed 
by our survey. The corresponding MF is reported in Table 7 
along with the value in Cha I for the same mass range. Both 
MF values are in very good agreement; the probability that 
they were drawn from the same distribution is 0.67. Bouy 
et al. (2006) surveyed low mass stars in the Upper Scorpius 
star-forming region using AO at the VLT. We consider only 
their subsample of 35 stars with spectral type M0-M6, cor- 
responding to —0.15-0.6 Mq at an age of 5 Myr according 
to the models of D'Antona & Mazzitelli (1997). The com- 
parison of this sample to ours (see Table 7) suggests a possi- 
ble possible higher incidence of binaries among the low-mass 
stars in Cha I compared to Upper Scorpius; the probability 
that the MFs in both regions, over this mass range, were drawn 
from the same distribution is 0. 1 1 . 

In a recent study of multiplicity over the separation range 
330-1650 AU, Kraus & Hillenbrand (2007) have found a sim- 
ilar trend of increasing binary fraction with increasing pri- 
mary mass in Upper Scorpius A, Taurus and Cha I. In partic- 
ular for Cha I, these authors determined binary fractions that 
increase from a few percents at 0.1-0.2M Q to —30% at 1.16- 
2.5 Mq. The observed dependence of binary fraction on pri- 
mary mass thus extends beyond the separation range probed 
by our survey. 

As can be seen from Figure 6, among our sample there are 
no triple systems in the lowest mass bin and there is a marginal 
increase in the fraction of triple systems with primary mass, 
although this trend is hampered by low number statistics in 
the highest mass bin. Such an increase is also observed by 
Correia et al. (2006) for their AO survey of wide binaries in 
various star-forming regions. 
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TABLE 7 

Multiplicity fractions for different regions and different masses 
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5.3. Mass ratio as a function of primary mass 

Figure 7 shows the mass ratio as a function of primary mass 
for all the multiple systems found. For this diagram, and 
the similar ones to come, each triple system contributes two 
points: one for the tight subsystem and one for the wide sys- 
tem. The mass ratios and separations for these are indicated 
in Table 4. Figure 7 clearly shows that as the primary mass 
decreases, the multiple systems are restricted to increasingly 
higher mass ratios. While such a trend could result from a 
selection effect, as lower mass ratio companions are more dif- 
ficult to detect around lower mass primaries, this is not the 
case here as our detection limits would have allowed us to de- 
tect systems with mass ratios much below the observed lower 
envelope at the low-mass end. The same tendency is also ap- 
parent among field multiple stars; this was recognized rather 
early based on a comparison of the mass ratio distributions 
of G-type stars, early-M stars, and late-M stars and brown 
dwarfs (e.g. Duquennoy & Mayor 1991; Fischer & Marcy 
1992; Close et al. 2003). 

Interestingly, the lower envelope of the mass ratio distribu- 
tion corresponds to a constant mass, which in fact roughly co- 
incides with the stellar/substellar boundary of ^0.075 Mq. In 
Cha I, according to our mass determinations, there are only a 
handful of companions lying just below the substellar bound- 
ary. However, given the uncertainties in mass and mass ratio 
determinations (see §3.2), it is possible that most or all of 
these companions have stellar masses. This low incidence of 
brown dwarf companions to stars is also found for G-M field 
stars, for which high-contrast imaging surveys revealed an in- 
cidence of brown dwarf companions of ^2-7% for the sep- 
aration range 25-1600 AU (Metchev 2006; Lafreniere et al. 
2007). This suggests that the process through which compan- 
ions form at separations greater than ~45 AU around stars 
more massive than ^0.1 M Q is inefficient in the substellar 
mass regime. However, the situation may be different for very 
low-mass stars (< 0. 1 Mq), for which several systems harbor- 
ing a brown dwarf companion are known. 7 

Figure 7 also reveals a dearth of near-equal mass binaries at 
the high mass end: the six systems with primary mass > 1 M Q 
have a mass ratio below 0.55. This is not the result of an ob- 
servational bias as high mass ratios would be easier to detect. 
The situation is similar for the high-mass binary systems in 

7 See the Very Low Mass Binaries Archive at 
http://paperclip.as.arizona.edu/ nsiegler/VLM_binaries/ 
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FIG. 7. — Mass ratios of the binary systems identified in our survey as 
a function of the primary mass. The solid line shows the mass ratio for a 
companion of 0.075 Mq. The dashed lines indicate the detection sensitivities 
reached for at least 50% of the targets for different angular separations. The 
dotted lines indicate a constant binding energy of 2 X 10 42 erg for various 
orbital separations, see §5.5 for more detail. Typical uncertainties on the 
estimated masses and mass ratios are ~15% and ~25%, respectively. 

Scorpius OB2, which also show a preference for low mass ra- 
tios (Kouwenhoven et al. 2005). The same trend is also found 
in the study of field G-type binary systems Duquennoy & 
Mayor (1991). The findings of Kraus & Hillenbrand (2007), 
for their combined Cha I, Taurus and Upper Scorpius sample, 
do not display this trend (see their figure 7); their sample con- 
tains ~12 binary systems with primary masses 1.16-2.5M Q 
and q > 0.75. However, of these ~T2 systems, more than half 
come from Upper Scorpius while only two come from Cha I. 
Moreover, one of these two Cha I systems (T26) clearly has 
a mass ratio much below 0.75: the 2MASS-based marginally 
resolved photometry used by Kraus & Hillenbrand (2007) to 
derive its mass ratio differs by ^3 mag from ours and from 
values reported elsewhere in the literature. Taking these fac- 
tors into consideration, their results are in agreement with 
ours, but may point toward a difference in high-mass binary 
properties between Cha I and Upper Scorpius. 

5.4. Mass ratio as a function of orbital separation 

Among our Cha I sample, two main features can be noted 
regarding the distribution of mass ratio versus separation. 
First, there is a clear lack of near-equal mass binaries (q > 0.5) 
at large separations (> 50 AU) even though such systems 
would be very easy to detect. Second, closer binaries tend 
to have mass ratios closer to unity. Both of these trends can 
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FIG. 8. — Mass ratio distributions for binaries with separation above (solid 
histogram) and below (dashed histogram) 50 AU. There is roughly an equal 
number of objects in both samples. 

be seen on Figure 8, which shows the mass ratio distributions 
for the subsets of small (< 50 AU) and large (> 50 AU) or- 
bital separations. The first mass ratio bin of the small separa- 
tion sample is the only one that could be significantly affected 
by incompleteness (see Table 2). By considering only sys- 
tems with a mass ratio over 0.2, where incompleteness should 
be unimportant, a Kolmogorov-Smirnov test indicates that 
the probabilities that the mass ratios of the large- and small- 
separations samples were drawn from a uniform distribution 
are 0.03 and 0.14, respectively, lending good statistical sig- 
nificance to the first trend mentioned above but only marginal 
significance to the second. The probability that the mass ratios 
for the small- and large-separation samples were drawn from 
the same distribution is only 0.02, again considering only sys- 
tems with a mass ratio above 0.2. 

Ratzka et al. (2005) looked at the flux ratio distributions 
for close (< 180 AU) and wide (> 180 AU) companions in 
p Oph (see their Fig. 13). As the flux ratio is a good proxy 
of mass ratio, their findings can be compared to ours. Be- 
fore proceeding further, we note that both features mentioned 
above are still be apparent if the separation between close and 
large separation systems is made at 180 AU instead of 50 AU, 
although the statistical significance is much reduced as there 
are only 12 systems with a separation greater than 180 AU. 
For the wide companions Ratzka et al. (2005) observed the 
same trend as us, namely, that wide pairs preferentially have 
large flux differences. However, the flux ratio distribution for 
their close pairs is more or less uniform, and does not show 
a peak near unity as for our sample. Similarly in Sco-Cen, 
Kohler et al. (2000) found a slight preference for wide pairs 
(> 190 AU) to have large flux differences, while the flux ra- 
tio distribution of close pairs is rather uniform. In Taurus, 
the wide pairs (> 180 AU) display a preference for large flux 
differences, and the close pairs also show a peak for flux ra- 
tios near unity (Kohler & Leinert 1998), very similar to our 
findings in Cha I. 

5.5. Orbital separation as a function of mass 

Figure 9 shows the projected separations of the binary sys- 
tems as a function of their total masses. It clearly shows a 
lack of wide systems at low masses; this is not an observa- 
tional bias as companions with larger separations are easier to 
detect. The observed decline of the maximum orbital separa- 
tion of Cha I binaries with decreasing mass is consistent with 
earlier findings for field systems (e.g. Burgasser et al. 2003; 
Close et al. 2003). The envelope of the separation distribu- 



FlG. 9. — Projected separations of the binary systems identified in our 
survey as a function of their total masses. The sizes of the symbols are pro- 
portional to the logarithm of the mass ratios of the binaries. The dashed line 
corresponds to a separation equal to 1500(M lot /MQ) 2 AU. Typical uncertain- 
ties on the estimated total masses are ~20%. 

tion of our Cha I sample is approximated reasonably well by 
1500(M tot /M Q ) 2 AU; only one binary, Hn 21W at 0.34 M Q 
and 880 AU, departs significantly from this trend. A very 
similar relation [1400(M tot /M Q ) 2 AU] was previously sug- 
gested for low-mass binary systems in the field (Burgasser 
et al. 2003), with an apparent break at ^0.6 Mq, above which 
wider pairs are found. For Cha I it is not possible, based 
on our data, to assess whether this relation is valid for to- 
tal masses greater than ^0.8 M Q since the corresponding or- 
bital separations are too close to or above our upper limit of 
~1000 AU. As the above relation between orbital separation 
and mass squared is reminiscent of a constant binding energy 
cutoff, we have calculated the binding energies of all binary 
systems. As suggested by the above trend, all systems (except 
Hn 21W) have binding energies above ^2 x 10 42 erg. As be- 
fore, this boundary is meaningful only for total masses below 
~0.8 Mq as our observations were much less sensitive to bi- 
naries with lower binding energies beyond this mass. Also, 
the uncertainty on this value is quite large, up to 50-60%, due 
to the uncertainties in masses, mass ratios, and projection ef- 
fects. 

The apparent lower boundary of binding energy could give 
rise to a trend such as the one noted earlier that less massive 
binaries preferentially have mass ratios closer to unity. In- 
deed, for low primary masses, the occurrence of companions 
of even lower masses could be suppressed as these would have 
binding energies that were too low. To investigate this possi- 
bility we have over plotted lines of constant binding energy 
(^2 x 10 42 erg) on figure 7 for four angular separations. Com- 
parison of these lines with the data indicates that for separa- 
tions below 1", 0.5", and 0.25", and primary masses above 
0.75 M Q , 0.3 Mq, and 0.2 Mq, respectively, systems with 
mass ratios much lower than the minima observed would still 
have binding energies above the minimum observed. Thus it 
does not appear that the envelope of the mass ratio-primary 
mass distribution is dictated by binding energy. Rather, this 
diagram seems to be better explained by a fixed minimum 
mass for companions. The binding energy appears to be re- 
sponsible only for the decrease of the orbital separation with 
decreasing mass. 

5.6. Stability of triple systems 

Tokovinin (2004), by studying a sample of triple systems 
from the Multiple Star Catalog, determined the following em- 
pirical stability criterion, P w /P t (l-e w ) 3 > 5, and an upper cut- 
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off on the periods ratios of P w /P t < 10 4 , where P w and P t are 
the periods of the wide and tight subsystems and e w is the ec- 
centricity of the wide subsystem. This criterion cannot be ver- 
ified for our triple systems as their true orbital separations and 
eccentricities are unknown. Nevertheless, the pseudo-periods 
calculated above (see Table 4) can provide some useful infor- 
mation. At first glance, all systems appear stable as they all 
have a pseudo-period ratio that satisfies the empirical stability 
criterion. Three systems (CHXR 28, CHXR 68A, and T26) 
have pseudo-period ratios well above the minimum value of 
5, indicating that they are almost certainly stable hierarchical 
triple systems. The other three systems, CHXR 9C, T31, and 
T39, have ratios of 12, 8.2, and 5.3, respectively, and thus 
could be unstable depending on their orbital orientations and 
eccentricities. 

We have performed a simple Monte Carlo simulation to in- 
vestigate the likelihood of stability of these three triple sys- 
tems. First we have randomly defined 10 5 fiducial stable 
triple systems filling uniformly the stability region of the e w 
vs log [P w /P t (l -e w ) 3 ] plane. The tight and wide pairs of each 
fiducial system were then randomly assigned an orbital orien- 
tation, phase, and eccentricity, and the corresponding projec- 
tion factors were calculated. The eccentricities of the wide 
pairs were drawn from a distribution following /(e) oc 2e. 
Then using the projection factors, a pseudo-period ratio was 
derived for each fiducial stable system. Of all fiducial sta- 
ble triple systems, only 6.3%, 4.1%, and 3% have a pseudo- 
period ratio smaller than 12, 8.2, and 5.3, respectively. Based 
on these fractions, it appears possible that one or more of 
CHXR 9C, T31, and T39 is indeed unstable; however we can- 
not reach a definite conclusion. 

5.7. Comparison with numerical simulations 

In this section we compare our observations mainly with 
the results of Goodwin et al. (2004a) and Delgado-Donate 
et al. (2004), which have done numerical simulations of 
the collapse and fragmentation of an ensemble of low-mass 
(~5 M Q ), dense, turbulent prestellar cores. Goodwin et al. 
(2004a) considered a low level of turbulence, and explored 
the effect of higher levels of turbulence in a subsequent paper 
(Goodwin et al. 2004b); they have followed the evolution of 
the system for 0.3 Myr. Delgado-Donate et al. (2004) consid- 
ered a level of turbulence about 20 times higher than Good- 
win et al. (2004a), and after a hydrodynamical phase lasting 
^0.5 Myr, they have followed the dynamical evolution of the 
systems for another 10 Myr using N-body simulations. Other 
differences between both sets of simulations include the ini- 
tial Jeans number, the initial density profile, and the equation 
of state. We caution the reader that, while these simulations 
are state-of-the-art and of very good quality, the limited range 
of core parameters they cover and the limited statistics they 
provide may impact the validity of our comparison to the en- 
tire population of Cha I, but some interesting observations can 
be drawn nonetheless. 

The simulations of Goodwin et al. (2004a) yield an overall 
MF of 0.26 and a CSF of 0.47. Apart from an unstable quin- 
tuple, they do not form systems of multiplicity higher than 4. 
In fact, the two quadruple systems they form would appear 
to us as triple or binary given our resolution limit. Good- 
win et al. (2004a) provide the separations and masses of all 
of their multiple systems, it is thus possible to interpret their 
results within the regime where our observations were sensi- 
tive. By treating systems with separations <16 AU as single 
and ignoring single stars less massive than 0. 1 M Q , we obtain 



a sample of 31 singles, 15 doubles, and 1 triple; this gives an 
MF of 0.34 and a CSF of 0.36. These do not differ signif- 
icantly from the quantities mentioned above for Cha I (0.27 
and 0.32 for the same separation limits). Among the high- 
order multiple systems in their sample, all but one of the 10 
binary subsystems have a separation less than 16 AU and half 
have a separation less than 10 AU. It will be interesting to 
see, with future radial velocity surveys or higher angular res- 
olution observations, whether many of the components of the 
multiple systems in Cha I are indeed unresolved binaries as in 
the Goodwin et al. (2004a) simulations. A higher level of tur- 
bulence leads to the formation of a larger number of stars per 
core (Goodwin et al. 2004b) and as a result, the average num- 
ber of components per multiple system is higher, although the 
MF is not affected much as a larger number of single stars 
are also ejected from the cores. The results of Goodwin et al. 
(2004b) indicate that the raw CSF may rise from 0.47 to 0.76 
when turbulence is increased by a factor of 5 compared to 
their baseline simulation of Goodwin et al. (2004a). However, 
applying the same procedure as above yields an MF of ^0.40 
and a CSF of ^0.42, which is still higher than the lower tur- 
bulence case but only by a small margin. Higher levels of tur- 
bulence, such as in the simulations of Delgado-Donate et al. 
(2004), may have an effect that can be readily compared with 
our results. In these simulations, even after 10.5 Myr, about 
1/3 of the multiple systems are quintuples or sextuples, which 
seems inconsistent with our observations even if a significant 
number of binary subsystems were unresolved given our reso- 
lution limit. Thus this excess of high-order multiples might be 
regarded as an indication that the level of turbulence in Cha I 
is well below that assumed by Delgado-Donate et al. (2004) 
in their simulations. 

An increase of the MF with primary mass, as observed for 
our Cha I sample (§5.2), is also found in the numerical sim- 
ulations of both Goodwin et al. (2004a) and Delgado-Donate 
et al. (2004); this trend was also recognized in the dynam- 
ical simulations of small clusters carried out by Sterzik & 
Durisen (2003). Furthermore, the increase of high-order mul- 
tiplicity with increasing primary mass we have noted earlier 
for our Cha I sample is also present in both models. To re- 
late the simulations to our observations, we have divided the 
systems formed in the simulations of Goodwin et al. (2004a) 
into the same mass bins as for our sample and the MF was 
calculated within each bin by treating binaries with separa- 
tion below 12 AU as single objects; the results are shown in 
Figure 6. The agreement with the observations is good, al- 
though two slight differences must be noted. First, the in- 
crease of the MF with mass is slightly steeper for the simula- 
tions than for the observations. Second, the MF in the lower 
mass bins are too low in the simulations. The overall behavior 
in the simulations of Delgado-Donate et al. (2004) is similar 
(see their Fig. 3), although the discrepancies just mentioned 
are more pronounced. As remarked by Delgado-Donate et al. 
(2004), the results of these simulations are dependent on the 
total mass of the cloud modeled. For instance, the simulation 
of a more massive core would lower the MF of 2-3 M Q stars, 
whereas a lower core mass could possibly raise the MF at the 
low-mass end. So perhaps a better agreement with observa- 
tions could be reached by simulating an ensemble of cores of 
different masses; this approach remains to be investigated. 

The simulations of Goodwin et al. (2004a) seem to repro- 
duce the trend we have observed in Cha I (Fig. 8) that large 
separation binaries preferentially have smaller mass ratios: 
all of their systems with separation larger than 30 AU have 
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q < 0.7 and all of their systems with separation below 30 AU 
have q > 0.4. In a similar vein, the simulations of Bate (2000) 
also predict that closer binaries should have more equal mass 
ratios, as observed in Cha I, owing to increased accretion on 
companions that form at smaller separations. 

6. CONCLUDING REMARKS 

We find that the young binary population in Chamaeleon I 
shares many trends with their field counterparts: the multi- 
plicity fraction declines and the mass ratios approach unity 
with decreasing primary mass, the mean and maximum sepa- 
ration of binary systems decrease with total mass, and brown 
dwarfs as companions to stars are rare. Of all the proper- 
ties we were able to investigate, the only clear difference be- 
tween the Cha I region and the field is the multiplicity frac- 
tion, which is lower in the field by a factor of at least ^ 1 .85 for 
the 16-1000 AU orbital range. Previous multiplicity surveys 
indicate that many properties of multiple systems, including 
their overabundance compared to the field, are similar in Cha I 
and in other dispersed associations, but more sensitive studies 
are needed to investigate possible differences among them. 
The findings to date seem consistent with the assertion that 
only a small fraction of the stars in the field comes from dis- 
persed star-forming regions, because if this were not the case 
a drop in multiplicity fraction with time without any effect on 
other properties of the multiple systems would be difficult to 
explain. However, these remarks may be somewhat premature 
given the relatively small number of multiples in our sample 
(36) and the lack of constraints on the very important orbital 
separation range below ~16 AU. In any case, it will be par- 



ticularly interesting to compare the same properties for multi- 
ple systems in high-density star-forming clusters, but this is a 
more difficult task given the larger distances to those clusters. 

Comparisons between multiplicity surveys and numerical 
simulations may also provide valuable insight into the star 
formation process. In particular, the fraction of high-order 
multiples could be related to the level of turbulence in the 
pre-stellar cores. The apparent paucity of such systems in our 
Cha I sample could be interpreted as evidence for a low level 
of turbulence in this region's cores. It may be that some of the 
binaries and triples in our sample harbor closer-in unresolved 
companions. Radial velocity surveys are needed to test that 
possibility. 
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APPENDIX 

EFFECT OF DIFFERENT EVOLUTION MODELS 

While it is not the goal of this work to compare the validity or the accuracy of the different evolution models available, it is 
worthwhile to investigate how a particular choice of model would influence the general trends found in our study. To address this 
issue, we have recalculated the masses of all primaries and companions using each of four different pre-main-sequence evolution 
models: D'Antona & Mazzitelli (1997), Siess et al. (2000), Palla & Stahler (1999), and Baraffe et al. (1998). When doing the 
analysis with each model, only the subset of objects falling within the effective temperature/magnitude range covered by that 
model was considered. The primary and companion masses were calculated in the way described in §3.2. For the models of 
Palla & Stahler (1999), magnitudes were computed from luminosities and effective temperatures using the method described in 
Appendix A of Testi et al. (1998); the tabulated magnitudes of the models were provided to us by F. Palla. 
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FIG. Al. — Multiplicity fraction as a function of primary mass, as cal- 
culated using different pre-main-sequence models. The dotted vertical lines 
mark the mass bins used. The masses and mass ratios were calculated using 
the models of D'Antona & Mazzitelli (1997) (diamonds), Siess et al. (2000) 
(squares), Palla & Stahler (1999) (circles), and Baraffe et al. (1998) (trian- 
gles); see §A for more detail. 
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FIG. A2. — Mass ratios as a function of projected separation, as calculated 
using different pre-main-sequence models. The masses and mass ratios were 
calculated using the models of D'Antona & Mazzitelli (1997) (plusses), Siess 
et al. (2000) (stars), Palla & Stahler (1999) (triangles), and Baraffe et al. 
(1998) (crosses); see text for more detail. 
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We concentrate on two diagnostic graphs: the MF as a function of mass, shown in Figure Al, and the binary mass ratio as a 
function of primary mass, shown in Figure A2. In computing the MF in the highest mass bin using the models of Baraffe et al. 
(1998), since these models do not cover the whole bin, we have assumed an upper cutoff of 6000 K. Figure Al shows that a very 
similar increase of the MF with mass is recovered using all models. The most important discrepancy occurs in the highest mass 
bin, where the point obtained using the models of Baraffe et al. (1998) is significantly below the other values. This lower MF 
value results from the higher ( 1 .5-2 x ) mass estimates obtained with the models of Baraffe et al. ( 1 998) for effective temperatures 
above ^3500 K; this shifts many single targets into the highest mass bin. Figure A2 indicates that a very similar distribution of 
mass ratio as a function of separation is also obtained using all models even though the primary masses can vary by as much 
as a factor of 2. Of the four sets of models considered, only those of D'Antona & Mazzitelli (1997) and Baraffe et al. (1998) 
cover the brown dwarf mass regime and a slight difference between them is worth mentioning. While the models D'Antona & 
Mazzitelli (1997) yield a handful of companions with masses slightly below the substellar boundary, the models of Baraffe et al. 
(1998), which typically yield higher primary masses, do not indicate any such companions. The lower envelope of the mass 
ratio distribution using the models of Baraffe et al. (1998) corresponds to a fixed mass of ~0.09 M Q rather than the 0.075 M Q 
value mentioned above for the models of D'Antona & Mazzitelli (1997). Despite the small differences noted above, our analysis 
indicates that the specific choice of evolution model has little effect on the general multiplicity properties of the sample. 

REFERENCES 



Ahmic, M., Jayawardhana, R., Brandeker, A., Scholz, A., van Kerkwijk, 

M. H., Delgado-Donate, E., & Froebrich, D. 2007, ApJ, 671, 2074 
Baraffe, I., Chabrier, G., Allard, E, & Hauschildt, P. H. 1998, A&A, 337, 403 
Bate, M. R. 2000, MNRAS, 314, 33 

Bate, M. R., Bonnell, I. A., & Bromm, V. 2003, MNRAS, 339, 577 
Bessell, M. S. 1991, AJ, 101, 662 

Bouy, H., Martin, E. L., Brandner, W., Zapatero-Osorio, M. R., Bejar, V. J. S., 

Schirmer, M., Huelamo, N., & Ghez, A. M. 2006, A&A, 451, 177 
Brandeker, A., Jayawardhana, R., Khavari, P., Haisch, Jr., K. E., & Mardones, 

D. 2006, ApJ, 652, 1572 
Burgasser, A. J., Kirkpatrick, J. D., Reid, I. N., Brown, M. E., Miskey, C. L., 

& Gizis, J. E. 2003, ApJ, 586, 512 
Burgasser, A. J., Reid, I. N., Siegler, N., Close, L., Allen, R, Lowrance, R, 

& Gizis, J. 2007, in Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & 

K. Keil, 427^141 

Chabrier, G., Baraffe, I., Allard, E, & Hauschildt, P. 2000, ApJ, 542, 464 
Close, L. M., Siegler, N., Freed, M., & Biller, B. 2003, ApJ, 587, 407 
Correia, S., Zinnecker, H., Ratzka, T., & Sterzik, M. F. 2006, A&A, 459, 909 
Cushing, M. C, Rayner, J. T., & Vacca, W. D. 2005, ApJ, 623, 1115 
Damjanov, I., Jayawardhana, R., Scholz, A., Ahmic, M., Nguyen, D. C, 

Brandeker, A., & van Kerkwijk, M. H. 2007, ApJ, 670, 1337 
D'Antona, F, & Mazzitelli, I. 1997, Memorie della Societa Astronomica 

Italiana, 68, 807 

Delfosse, X., Beuzit, J.-L., Marchal, L., Bonfils, X., Perrier, C, Segransan, 
D., Udry, S., Mayor, M., & Forveille, T. 2004, in Astronomical Society 
of the Pacific Conference Series, Vol. 318, Spectroscopically and Spatially 
Resolving the Components of the Close Binary Stars, ed. R. W. Hilditch, 
H. Hensberge, & K. Pavlovski, 166-174 

Delgado-Donate, E. J., Clarke, C. J., Bate, M. R., & Hodgkin, S. T. 2004, 
MNRAS, 351, 617 

Duchene, G., Bouvier, J., & Simon, T. 1999, A&A, 343, 831 

Duquennoy, A., & Mayor, M. 1991, A&A, 248, 485 

Fischer, D. A., & Marcy, G. W. 1992, ApJ, 396, 178 

Fitzpatrick, E. L. 1999, PASP, 111, 63 

Frink, S., Roeser, S., Alcala, J. M., Covino, E., & Brandner, W. 1998, A&A, 
338, 442 

Ghez, A. M., McCarthy, D. W., Patience, J. L., & Beck, T. L. 1997, ApJ, 481, 
378 

Ghez, A. M., Neugebauer, G., & Matthews, K. 1993, AJ, 106, 2005 
Goodwin, S. P., Whitworth, A. P., & Ward-Thompson, D. 2004a, A&A, 414, 
633 

— . 2004b, A&A, 423, 169 

H0g, E., Fabricius, C, Makarov, V. V., Urban, S., Corbin, T., Wycoff, G., 

Bastian, U., Schwekendiek, P., & Wicenec, A. 2000, A&A, 355, L27 
Jayawardhana, R., & Ivanov, V. D. 2006, Science, 313, 1279 
Kenyon, S. J., & Hartmann, L. 1995, ApJS, 101, 117 
Kohler, R., Kunkel, M., Leinert, C, & Zinnecker, H. 2000, A&A, 356, 541 
Kohler, R., & Leinert, C. 1998, A&A, 331, 977 

Kohler, R., Petr-Gotzens, M. G, McCaughrean, M. J., Bouvier, J., Duchene, 
G, Quirrenbach, A., & Zinnecker, H. 2006, A&A, 458, 461 



Konopacky, Q. M., Ghez, A. M., Rice, E. L., & Duchene, G. 2007, ApJ, 663, 
394 

Kouwenhoven, M. B. N., Brown, A. G. A., Zinnecker, H., Kaper, L., & 

Portegies Zwart, S. F. 2005, A&A, 430, 137 
Kraus, A. L., & Hillenbrand, L. A. 2007, ApJ, 662, 413 
Kraus, A. L., White, R. J., & Hillenbrand, L. A. 2005, ApJ, 633, 452 
— . 2006, ApJ, 649, 306 
Kroupa, P. 1995, MNRAS, 277, 1522 
— . 1998, MNRAS, 298, 231 
Lada, C. J. 2006, ApJ, 640, L63 

Lafreniere, D., Doyon, R., Marois, C, Nadeau, D., Oppenheimer, B. R., 
Roche, P. F, Rigaut, F, Graham, J. R., Jayawardhana, R., Johnstone, D., 
Kalas, P. G., Macintosh, B., & Racine, R. 2007, ApJ, 670, 1367 

Leinert, C, Zinnecker, H, Weitzel, N., Christou, J., Ridgway, S. T., Jameson, 
R., Haas, M., & Lenzen, R. 1993, A&A, 278, 129 

Lodieu, N., Hambly, N. C, Jameson, R. E, & Hodgkin, S. T. 2008, MNRAS, 

383 1385 
Luhman, K. L. 2004, ApJ, 602, 816 

Luhman, K. L., Wilson, J. C, Brandner, W., Skrutskie, M. F, Nelson, M. J., 
Smith, J. D., Peterson, D. E., Cushing, M. C, & Young, E. 2006, ApJ, 649, 
894 

Masciadri, E., Brandner, W., Bouy, H., Lenzen, R., Lagrange, A. M., & 

Lacombe, F. 2003, A&A, 411, 157 
Mathieu, R. D., Ghez, A. M., Jensen, E. L. N., & Simon, M. 2000, Protostars 

and Planets IV, 703 
Metchev, S. A. 2006, PhD thesis, California Institute of Technology, United 

States - California 
Padgett, D. L., Strom, S. E., & Ghez, A. 1997, ApJ, 477, 705 
Palla, E, & Stahler, S. W. 1999, ApJ, 525, 772 

Petr, M. G., Coude Du Foresto, V, Beckwith, S. V. W., Richichi, A., & 

McCaughrean, M. J. 1998, ApJ, 500, 825 
Ratzka, T, Kohler, R., & Leinert, C. 2005, A&A, 437, 61 1 
Reid, I. N., & Gizis, J. E. 1997, AJ, 113, 2246 

Rousselot, P., Lidman, C, Cuby, J.-G., Moreels, G., & Monnet, G. 2000, 

A&A, 354, 1134 
Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 
Simon, M., Close, L. M., & Beck, T. L. 1999, AJ, 117, 1375 
Skiff, B. A. 2007, VizieR Online Data Catalog, 1, 2023 
Sterzik, M. E, & Durisen, R. H. 1998, A&A, 339, 95 
Sterzik, M. E, & Durisen, R. H. 2003, A&A, 400, 1031 
Testi, L., Palla, F, & Natta, A. 1998, A&AS, 133, 81 

Tokovinin, A. 2004, in Revista Mexicana de Astronomia y Astrofisica 

Conference Series, Vol. 2 1 , Revista Mexicana de Astronomia y Astrofisica 

Conference Series, ed. C. Allen & C. Scarfe, 7-14 
Whittet, D. C. B., Prusti, T, Franco, G. A. P., Gerakines, P. A., Kilkenny, D., 

Larson, K. A., & Wesselius, P. R. 1997, A&A, 327, 1 194 
Zacharias, N., Urban, S. E., Zacharias, M. I., Wycoff, G. L., Hall, D. M., 

Monet, D. G., & Rafferty, T. J. 2004, AJ, 127, 3043 



